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t^j- . Abstract 

CO; 

We consider the effects of metastable charged sparticles on Big-Bang Nucleosynthesis (BBN), in- 
cluding bound-state reaction rates and chemical effects. We make a new analysis of the bound 
states of negatively-charged massive particles with the light nuclei most prominent in BBN, and 
present a new code to track their abundances, paying particular attention to that of 7 Li. Assum- 
ing, as an example, that the gravitino is the lightest super symmetric particle (LSP), and that the 
^ ■ lighter stau slepton, f\, is the metastable next-to-lightest sparticle within the constrained minimal 

supersymmetric extension of the Standard Model (CMSSM), we analyze the possible effects on the 
standard BBN abundances of f\ bound states and decays for representative values of the gravitino 
mass. Taking into account the constraint on the CMSSM parameter space imposed by the discovery 
of the Higgs boson at the LHC, we delineate regions in which the fit to the measured light-element 
abundances is as good as in standard BBN. We also identify regions of the CMSSM parameter 
space in which the bound state properties, chemistry and decays of metastable charged sparticles 
can solve the cosmological Li problem. 
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1 Introduction 



The agreement of standard Big-Bang Nucleosynthesis (BBN) calculations with the measured 
abundances of the light elements imposes important constraints on scenarios for new physics 
that involve massive metastable particles [1] - [17]. If these particles are neutral, only the 
effects of particles produced in the showers following their decays need to be taken into 
account, but in the case of negatively-charged metastable particles X~, the formation of 
(AX') bound states should also be considered, and are very important [23ll26l429ll3 1 U34f436l 

susu. 

The emergence of the primordial 'lithium problem' adds cosmo logical motivation to the 
studying the effects of metastable particles on BBN. As reviewed in refs. [39 | l4*4" t l4*8 | l4"9]. 
WMAP and other observations [50] have determined precisely the cosmic baryon density 
and thus pinned down the one free parameter of standard BBN |51j. Using this as an 
input, BBN makes precise predictions for light-element abundances, and those of deuterium 
and 4 He are in good agreement with obervations. But the BBN expectations for 7 Li/H 
based on the WMAP baryon density are higher than the observed abundances by factors 
of 2 — 4, amounting formally to a 4 — 5cr discrepancy; this is the cosmological lithium 
problem J52j[53]. Nuclear uncertainties [53H56] and/or resonances are all but excluded as 
solutions to the problem [5*6TI6"T] . There is the possibility that depletion plays a role in altering 
the 7 Li abundance [Q2\. However, these solutions typically have difficulty in explaining the 
thinness of the 7 Li plateau [63J as well as the observation of 6 Li in some halo stars [61]. The 
temperature scale used in the 7 Li abundance determination has also been considered [65|l6"6] 
and it seems unlikely that a significant change in the 7 Li abundance is possible within 
reasonable uncertainties in the effective temperature. We note, however, recent observations 
of lithium in the interstellar medium of the metal-poor Small Magellanic Cloud test these 
systematics and are consistent with the halo-star results |67] . Thus, the cosmological lithium 
problem seems increasingly likely to be real, and to point to new physics during or after BBN. 

In a previous paper |43] . we extended analyses of the effects of particle showers in the 
decays of metastable particles to include the most relevant uncertainties in nuclear reaction 
rates. We applied our analysis to scenarios within the constrained minimal supersymmetric 
extension of the Standard Model (CMSSM, see Appendix B for its specification) in which the 
lightest neutralino x is the lightest supersymmetric particle (LSP), and the heavier, neutral 
gravitino is metastable. Not only did we find regions of this CMSSM parameter space where 
the cosmological light-element abundances agreed with the measured values at least as well 
as in standard BBN, but we also identified regions of this CMSSM parameter space where 
the cosmological 7 Li problem is alleviated and even potentially solved. In this paper, we 
extend the analysis of [13] to include the (AX~) bound-state effects expected in the case of 
a negatively-charged metastable particle X~ . 

Bound-state effects were also discussed in [26], and here we update and supersede that 
analysis incorporating qualitatively and quantitatively new rates and processes that were not 
available at the time. To this end, we first review our calculations of bound-state properties, 
and then turn to their effects on BBN. These include calculations of (1) bound-state recombi- 
nation, which fixes the abundances of various exotic 'ions' such as (pX~), ( 4 HeX~), ( 7 BeX~), 
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etc., and (2) bound-state catalysis, which causes additional changes in light-element produc- 
tion and destruction rates beyond the non-thermal reactions considered in jl3]. Our new 
calculations of bound-state properties such as binding energies and charge radii are in rea- 
sonable agreement with other work, and we use them to discuss the effects of uncertainties 
in the nuclear inputs. For this purpose, we have compiled a complete and up-to-date list 
of the relevant reactions, tabulated below. We have verified, using a simple driver code, 
that our recombination rates give (pX~) and ( 4 HeX~) abundances in good agreement with 
previous results [23]. We have then updated the BBN code used in [261S01S3] to include the 
recombination and catalysis rates, in a more complete, accurate and systematic way than 
previously. 

As an example of the application of our code, we consider the case of a supersymmetric 
model in which the gravitino is the LSP, and the lighter stau slepton, fx is the metastable 
NLSP. We work within the framework of the CMSSM, and seek regions of its parameter 
space where the consistency with the measured light-element abundances of standard BBN 
calculations is at least maintained, and also look for regions where the cosmological 7 Li 
problem may be alleviated or even solved. We find that this is possible for generic values of 
the CMSSM parameters where the lifetime of the NLSP ~ 10 3 s (as in p~5l[16]), and that 
there are more extended regions of parameter space where the cosmological 7 Li is at least 
no worse than in standard BBN calculations. 

The paper is organized as follows. In Section 2 we discuss the relevant properties of 
X~ bound states, including the Coulomb radii of several nuclides and our three-body model 
for the (a + a + X) system, and various choices for binding energies. In Section 3, we 
discuss relevant nuclear interaction rates involving bound states, and our implementation of 
them in the BBN network. In Section 4, we briefly describe the chemical reactions involving 
bound states. In Section 5, we introduce the supersymmetric framework we use to produce 
our numerical results. In particular, we consider models where the lighter stau slepton 
is our candidate for the metastable charged particle X~, and results for stau lifetimes are 
summarized in Section 6. Our main results are given in Section 7, which includes a discussion 
of the light element abundance observations and the abundances we find from BBN with 
stau bound states. Our results are summarized in Section 8. 

2 Bound-State Properties 

Before considering the impact of a new electromagnetically-charged particle on BBN pre- 
dictions, we first study the properties of its bound states with light nuclei. To this end, we 
solve the time-independent Schrodinger equation with an interaction potential given by the 
Coulomb potential between a finite-sized nuclide and a point-like X~ particle. The Coulomb 
potential is determined by the charge and the rms charge radius of the nuclide of interest. 
Here we adopt the latest rms charge radii measurements of He, Li and Be isotopes, and 
assume the charge distribution of the nuclide to be Gaussian 0. 

1 As long as one reproduces the rms charge radius, the detailed form of the charge distribution is not 
important for our purposes. 
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When solving the Schrodinger equation, we first define dimensionless variables for the 
energy and a typical length scale. This makes the equations a function of a single parameter, 
the ratio of the rms charge radius R c to the Bohr radius Rb- We have verified that our 
solutions interpolate smoothly between Coulombic bound-state energies at small R c /Rb and 
harmonic oscillator energy levels at large R c /Rb, and that our numerical solutions match 
analytic solutions. 

Numerical results for bound states for the nuclides of interest for BBN are shown in 
Table [TJ We note that the rms charge radii of many of the nuclides considered have been 
determined experimentally. However, in other cases we must rely on phenomenological 
estimates, and Table Q] gives some ranges in these cases. In particular, the charge radii 
for the two nuclei relevant for crossing the A = 8 divide, namely 8 Be and the first excited 
state of 9 Be, are not known experimentally. Given that 8 Be is a barely bound state of two a 
particles ( 4 He nuclei) that are hardly touching, it is expected that the rms charge radius of 
8 Be should be close to twice the rms charge radius of 4 He, i.e., i? Cj8 = 3.362 fm. This estimate 
is in good agreement with the value given in [75], namely R c ^ = 3.39 fm, after correcting for 
the poor binding energy determination in [77]. For comparison, we also consider the value 
R c .s = 2.50 fm [76]; an estimate chosen close to the 9 Be charge radius, though likely an 
underestimate. 

The ( 8 BeX~) state can play an important role in primordial 9 Be synthesis, if the state 
exists and has appreciable abundance. As we will see in the next section, the significance 
of ( 8 BeX~) depends not only to the qualitative issue of whether this system is bound, but 
also on the value of the binding energy. For this species, therefore, we have gone beyond the 
2-body model that has been used to calculate the other binding energies appearing in Table 
[2j To analyze the bound state properties of the 3-body system (a + a + X), we have utilized 
the Variational Monte Carlo (VMC) method [78] . This method assumes a suitable form 
of wavefunction with parameters to be determined by Monte Carlo variation. Parameters 
are randomly chosen after some initial guess; the expectation value of the Hamiltonian is 
then computed and, if it is lower than the initial expectation value, the current parameters 
are set as the adopted values. This variation is repeated until some convergence criteria 
are met, such as small (or no) changes in the expectation value of the Hamiltonian with 
decreasing step size. The complexity of the adopted wave function can be increased, with 
added parameters, allowing for another test of convergence. Note that, as with the 2-body 
case, the 3-body binding energy is expressed with respect to free 8 Be+X; this is larger than 
the binding relative to free a + a + X by m( 8 ' free Be) — 2m a = —Qs,iree = 92 keV. 

We have adopted wavefunctions with the forms of exponentials with Gaussian cutoffs 
and extra Gaussian terms with decreasing dispersion with respect to the cutoff scale. The 
method was validated with the three-body systems of the neutral He-atom and the negatively 
charged H-atom, using finite-sized Coulomb potentials with the p and a rms charge radii. 
The VMC method reproduced the binding energies for these Coulomb-only systems quite 
quickly, agreeing with the observed values. For the a + a + X system, we similarly adopt 
finite-sized Coulomb interactions, plus an added nuclear a + a interaction from [79]I80]. The 
adopted nuclear a + a potential reproduces elastic scattering data and the 8 Be ground state 
resonance energy. We find that 8 Be + X is bound, but with B$ = 492 ± 50 keV. This is 



3 



Table 1: Properties of X~ bound states with relevant nuclides. The Table lists relevant 
nuclides together with their (unbound) masses and rms charge radii R c in Coulombic pa- 
rameterizations of the potentials. Bound-state binding energies Ba come from our 2-body 
calculations based on the given charge radii (except where otherwise noted). Experimental 
values and ranges of R c are listed, where available for some nuclides, and ranges of theoretical 
estimates of R c for other nuclides. 



Nuclide 


Mass (amu) 


Rc (fm) 


B A (MeV) 


!H 


1.00782503 


0.8750 ± 0.0068 [68| 


0.02493 


2 H 


2.01410178 


2.1303 ±0.0010 [69] 


0.04879 


3 H 


3.01604927 


1.63 ±0.03 [70] 


0.07264 


3 He 


3.01602931 


1.9506 ±0.0014 [71] 


0.2677 


4 He 


4.00260325 


1.681 ± 0.004 [7J 


0.3474 


6 Li 


6.01512228 


2.517 ±0.030 [73] 


0.8000 


7 Li 


7.01600405 


2.39 ±0.030 [73] 


0.8893 


7 Be 


7.01692925 


2.647 ±0.015 [71] 


1.2879 


8 Be 


8.00530509 


3.390 [75] 


1.1679 






2.50 [76] 


1.408 [76] 






N/A* 


0.492 


9 Be 


9.01218213 


2.519 ±0.012 [7J 


1.4699 


9 Be* 


9.01398998 


2.519 


1.4700 






2.880 [77] 


1.3527 






3.390 


1.2173 


8 B 


8.037675026 


2.65 


1.8547 



*Our result for Ba in this case is based on our three-body calculation. 

more fragile than found when assuming a two-body ( 8 BeX~) structure with an appropriate 
rms charge radius. The Coulomb repulsion of the two a particles loosens the binding of the 
three-body system. 

For ( 8 BeX~) to allow substantial 9 Be production, two conditions must be satisfied: 

1. Production of ( 8 BeX~) must be possible and effective, and 

2. Production of 9 Be must proceed resonantly through the first excited state in ( 9 Be*X~), 
that is, ( 8 BeX-) +n -»■ ( 9 Be*X~) -> 9 Be + X~ [76]. 

The requirement that ( 8 BeX~) production is possible merely demands that this state is 
stable, as we and others have found. But, as we will see, in the early Universe ( 4 HeX~) is 
the dominant bound state. Consequently, ( 8 BeX~) production occurs via ( 4 HeX~) ± 4 He — > 
( 8 BeX~) + 7. This channel is only effective if it has Q > 0, which demands that 

B 8 > Bf» = B A - Q 8 , frce = 0.439 MeV . (1) 

We find that all estimates of B$ satisfy this constraint, though our 3-body result does so by 
a much smaller margin than the others. 
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The second requirement for 9 Be production depends on the position of the first excited 
state of ( 9 Be*X~), to which we now turn. In ordinary 9 Be, this is a cluster state that is 
poorly described by shell-model calculations. Its structure is that of two a particles and a 
neutron and, to first order, this state is the 8 Be + n ground state. Therefore, a first guess for 
the rms charge radius would be R Cj g* = -R c ,s- However, the presence of the neutron impacts 
the structure, and the excited state of 9 Be should have a radius larger than the ground state. 
This constrains our estimate of the impact of the neutron on the rms 9 Be charge radius to 
the range R c ,9* £ (2.519,3.39) fm, where the low value is just the ground-state value R c ,9i 
and the upper value is the estimate of R c $ given in [75]. This range includes the result for 
the rms charge radius given in [77], namely -R c ,9* = 2.88 fm. If one assumes the same relative 
shift in the rms charge radius as for 8 Be, one finds R c $* = 3.11 fm. However, it may be 
larger, given that the level energy calculated is less accurate in the 9 Be case than in the 8 Be 
case. 

9Be* 

8Be+n fflBriO+n 8Be+n+X 

0.110 a+q+n ^otseAj+n q+q+n+X - 



^MXXT (4 HeX) + qT na00Q- 0092 

u.uuu -0.400 [ 



-0.347 



Kamimura et al _-j Qyg (9B6*X) 
Qg0 Pospelov - 



-1.316-1.360-1.107 

-1.574 



1 .470 MeV 



(9BeX) 
-3.044 



Figure 1: Level schemes calculated for unbound and bound X~ states in the 9 Be system. 
State labels appear above the corresponding level. Level values shown are based on Table 1, 
and all levels are shifted such that zero energy corresponds to the free-particle state. In the 
case of ( 8 BeX~) , we show the result of our 3-body calculation (a blue band whose the width 
corresponds to the uncertainty), as well as the 2-body calculations of Kamimura et al. J73]/ 
and of Pospelov |76|/ (labelled). In the case o/( 9 Be*X~), we consider a range of level values 
including our and Pospelov's estimates. 

Fig. [TJ shows the level structure for the 9 Be system when either in the ordinary unbound 
state (left, [81J ) or bound to an X~ particle (right, our calculations). The zero-point of 
energy is taken to be that of unbound, free particles. Thus, for example, unbound 8 Be + n 
lies at — Qs.free = +0.092 MeV. The level positions for each of the bound states are shifted 
relative to the corresponding unbound case due to the binding energy: level i lies at energy 
E(ix) = Eij ree — Bi. Thus, level spacings are shifted by differences in binding: E^ x \—Eu X ) = 
Ejfi ee — -E^free — Bj + Bi. For example, the ( 9 Be*X) excited state lies above the ground state 
by an amount that is the sum of the unbound 9 Be level spacing 1.684 MeV, minus the 
difference B 9 * — B 9 . If the bindings are the same, then the first excited state spacing is also 
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the same as in the ordinary case. But if the excited state is more weakly bound, then the 
level spacing is larger than in ordinary 9 Be. 

The three strips for the ( 8 BeX~) bound state represent the three values as determined by 
this work (highest strip), as well that of [751176] (the lower two strips, as labelled). The thick 
box for the excited ( 9 Be*X~) state corresponds to the range of possible charge radii -R Cj 9* G 
(2.519, 3.39) fm. Even with the broad range of possible rms charge radii for the first excited 
state in 9 Be, one can see that the excited state probably lies substantially below the ( 8 BeX~) 
entrance level. If this is so, this drives the reaction away from resonance, substantially 
reducing or even eliminating the mechanism discussed in [76], and thus suppressing 9 Be 
production via bound states. Ultimately, however, a true determination of the first exited 
state position would require a four-body calculation of (aanX~). This challenging work 
has not been done, and lies beyond the scope of our paper. Note also that our 2-body and 
3-body calculations of ( 8 BeX~) gave significantly different binding energies, underscoring 
the importance of detailed calculations for these states with relatives high Z. Consequently, 
we cannot exclude the possibility that the ( 8 BeX~) entrance channel will be near resonance. 

Clearly, the situation regarding 9 Be is uncertain, reflecting the poor state of knowledge 
regarding the ( 8 BeX~) and ( 9 BeX~) nuclear properties. In this paper, our approach is to 
illustrate the ability of 9 Be to constrain supersymmetric models, within the most optimistic 
scenario in which the resonant production discussed in [76] occurs. Thus, for most of our 
calculations we adopt this (constant) resonant rate for ( 8 BeX~) + n — > 9 Be + X~ [76J. 
However, the reader should bear in mind that the resulting 9 Be abundances and resulting 
constraint therefore represent a 'most optimistic' scenario. Consequently, we also make 
comparisons with calculations in which this production channel is suppressed, considering 
the cases: 

1. Bound state structure and resonant rate from |76j. 

2. ( 8 BeX~) is more tightly and more weakly bound. 

3. The resonant rate set to zero for ( 8 BeX~) 9 Be + X~ . 

It goes without saying that there is a pressing need for precise and accurate calculations of 
( 8 BeX~) and ( 9 BeX~) properties. 

3 Bound-State Reactions and Abundances 
3.1 Formalism 

In looking at the effects of bound states of X~, we must track the abundances of its bound 
states with various nuclei, e.g., (pX~), ( 4 HeX~) and ( 7 BeX _ ). For this purpose, we need to 
incorporate reactions that affect these bound-state abundances, namely: (1) recombination 
and photodissociation processes on each nuclear species i, such as 

i + X- <-> (zX")+7 , (2) 
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and (2) charge-exchange processes between (iX ) bound states and other nuclides j, such 



as 



(iX-)+j<*i+(jX-) . (3) 

We refer to these processes collectively as 'bound-state chemistry', and solve the rate equa- 
tions for the processes (J2]) and ([3]) to determine the corresponding chemical abundances. 

Denoting the total X~ number density by nx- = Yx-Ub, we decompose it into the free 
and bound abundances, respectively %- ifree and n^ x -), with j e p,d,. . .. To remove the 
effect of cosmic expansion, as usual we follow the evolution of the 'mole fractions' 



n, 



n B 



(4) 



where the states i include: (a) ordinary unbound, free nuclei, (b) X~ bound states, and (c) 
free X~ . Thus, we treat these in a manner completely parallel to the usual BBN accounting 
for ordinary (unbound) nuclides. Note that the total abundance of a nuclear species i sums 



i,tot 



^i.free ~t~ Y( 



(iX~), 



while the total abundance of X is 



its unbound and bound states Y it 

Yx~,tot — Yx-,fvee + Y(iX~)- 

For two-to-two reactions of the form ah — > cd, the reaction rate per unit volume is 
n a nb (cab^cdv) , with ((Tab^cdv) the appropriate thermally-averaged rate coefficient. The re- 
action rate per particle a is thus 



Tabbed — {o~ab— > cdV) U b = N Am ((Tab^cdV) pBY b = A a 

b^rcd 



(5) 



where A^Avo = l/ m « is Avogadro's number, and \ a b-*cd — A^Avo (<Jab->cdv) is the form in which 
thermonuclear rates are normally tabulated. Thus the total rate per target b nucleus is 

^ab\cd = ^ab-^cd + T ' (aX)b->cdX = (o ab^cdV) ^B^free + (o~(aX)b^cdXv) /iB^aX) , (6) 

which can be substantially larger than in the ordinary case if the catalyzed rate coefficient 
has a large enhancement and if there is a substantial (aX) abundance. 

With these definitions, the evolution of bound state (iX), with i 6 p,d, . . . can be ex- 
pressed as a sum over several kinds of processes: 



dt 



Y 



(tX) 



d_ 
dt 



Yt 



(iX) 



chcm 



(7) 



decay 



The bound-state chemistry reactions do not change the type of nuclides in the initial state, 
and are 



d_ 

Of 



Y (i 



(iX) 



chcm 



^ Y( y i X )+T i X ^ 1 (i X )Y x ^ cc + s y^ j r 



i(jx)^j(ix)Y(jx) 



Bound-state nuclear reactions have final-state nuclides different than those in the initial 
state, and take the form 



d_ 

dt 



Y 



(iX) 



Y v i(iX)^ktx Y {lX) + r 



e(kx)^j(ix)Y(kx) 



(9) 



kf. 



kf 
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where the last term includes only those reactions that produce bound (iX ) rather than free 
i. Finally, the decays of X~ with lifetime r% destroy bound states: 



d~t Y{lX) 



(10) 



decay 



where the decay rate T x = 1/tx- 
Turning to free X~ , we have 



d_ 
dt 



X,free 



7» Y W 



chem 



7» Y W 



-T X Y: 



X J-X, free 



The bound-state contribution to the evolution of a species of unbound, free nuclei i is given 
by 



d 



BS 



d 

~ dt Y[ ' x) 



chem 



d 



decay 



2^ ^i{jX)^k£X J ^j.free + T 
v j J U 



e(kx)^ijxY{kx) , 



'12) 



where the last term includes only bound-state reactions that produce free i in the final state. 



3.2 Reaction Rates 

Tables [2] and [3] summarize the treatments of bound-state chemistry and nuclear rates in our 
work and in the recent literature. Entries for recent literature are given to the best of our 
knowledge; in some cases full details were not given in the published papers. Blank entries 
mean that to the best of our knowledge no rate was assigned, effectively setting the rate to 
zero. 

Wherever possible, we adopted the most up-to-date chemistry and nuclear rates from the 
literature. In many cases, rates were not available for channels we wished to examine. Thus 
we adopted simple rules to estimate the needed rates from published ones; these cases are 
identified in the Tables. In the case of bound-state chemistry, we adopted recombination 
rates for nuclide i using the scaling cT rec oc ZfB( iX )- 

As we see below, bound state chemistry strongly favors ( 4 HeX~) production, which 
essentially locks up all of the X~, for the cases of physically interesting abundances where 
Yx- < Y a . Consequently, nuclear reactions involving ( 4 HeX~) are the most important. On 
the other hand, (pX~) reactions are relatively unimportant due to the small abundance of 
this state, and (dX~) and (tX~) have negligible effects. 

In general, bound states enhance nuclear rates. This is in part because they reduce the 
Coulomb barrier, to which the rates are exponentially sensitive. For nuclides and channels for 
which bound-state nuclear rates were not available in the literature, we estimated the rates 
assuming this is the only source of perturbation. In these cases, we adopted the ordinary 
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thermonuclear rates, but with a Gamow penetration factor appropriate for a nucleus of 
effective charge Zf^ x , = Z { — 1. 

In many cases, bound states also enhance nuclear channels through catalysis effects, 
which may be described as follows [23]. Consider the important example of catalyzed 6 Li 
formation, 

( 4 HeX") + 6 Li + X" . (13) 

The corresponding ordinary process is the 4 He + d — > 6 Li + 7 radiative capture reaction, 
which is suppressed because it must proceed through the E2 mode. The bound-state rate 
does not require the emission of a photon and is substantially larger than the ordinary rate 
in the typical situation in which the ( 4 HeX~) abundance is large. 

We have included rates for the formation and processing of ( 8 BeX~) which, as we will see, 
can lead to substantial 9 Be production under optimistic nuclear physics assumptions. We 
also have included rates for 10 B and n B production via rates involving (aX~). We find that 
boron production is indeed increased over the (very small) standard BBN level. However, 
the B/H abundance always remains many orders of magnitude below the levels seen in halo 
stars. Thus we find that boron is not a promising signature of decaying particle effects. 

We have also studied whether reionization by the emitted XLya photons could inhibit 
the net rates for the NLSP recombination reactions A + X~ — > (AX~) + 7, as is the case in 
ordinary hydrogen recombination. As discussed in Appendix A, we find that, whereas the 
optical depth for reionization by XLya photons emitted by NLSP recombination is much 
smaller than that for ordinary hydrogen and helium recombination, it is still very large, so 
that the net rate of recombination might be very suppressed. However, as we also show in 
Appendix A, XLya: photons Compton scatter rapidly off free electrons. This rapidly degrades 
the energies off resonance, so that they are ineffective for reionization. We conclude that 
NLSP recombination to the ground state proceeds unimpeded, unlike the case of ordinary 
hydrogen and helium recombination. 

4 Bound-State Chemical Effects 

We present later results from a code that treats self-consistently the bound states as separate 
nuclei, which then can have their own set of bound-state chemical and nuclear reactions with 
other species. As a warm-up exercise, we first present some results with catalysis effects 
turned off, and so only incorporate bound-state chemistry, i.e., recombination onto bound 
states. We include decays as part of the chemistry, i.e., decays remove free and bound X~ , 
but we turn off nonthermal decay effects. This exercise tests our code and illustrates the 
interplay between recombination and charge transfer. 

For this purpose, we choose an initial X~ abundance Y^ tot = 10~ 2 , which is typical for 
interesting supersymmetric models, and we vary the lifetime Tx, to show the sensitivity to 
this parameter. 

In Figs. |5] - IH we show the abundances Yi = rti/n-Q for both bound and free species, 
as functions of the temperature T. The solid black line corresponds to the abundance of 
free X~'s, whereas the other solid lines are the abundances of the bound states, as labeled 
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Table 2: Summaries of the treatments of bound-state chemistry rates assumed in our work 
and in the recent literature (I) 











CEFOS 


Bailly 


Pospelov 


Kamimura 


This 




Reaction 




2006 


et al 2009 [37] 


et al 2008 [3J] 


et al 2009 [75] 


Work 
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1 


simple scaling 


Kamimura 09 [75] 


estimated 




Pospelov [36] 


c?7 


x~ 


- -> (dX-) + 


7 


simple scaling 








scaled Pospelov [36] 


t + 


x~ 


->(**-) + 


7 


simple scaling 








scaled Pospelov [36] 


3 He4 


x~ 


~ -> ( 3 HeX-; 


+ 7 


simple scaling 








scaled Pospelov [36] 


a 7 


x~ 


- -> (aX-) 4 


7 


simple scaling 




Pospelov 07 [53] 




Pospelov [35] 


6 Li7 


x~ 


-> ( 6 LiX-) 


+ 7 


simple scaling 








cr rcc oc Z 2 B scaling 


7 Li7 


x 


- -> ( 7 LiX-) 


+ 7 


simple scaling 








er rcc oc Z 2 B scaling 


7 Bc7 


x 


- -> ( 7 BeX") 


+ 7 


simple scaling 




Bird 08 [28] 




Bird 08 [28] 


{ 8 BcX-) 


+ 7 -> 8 Be + 


X~ 


simple scaling 








cr r cc oc Z 2 B scaling 


9 Be + 


x 


- -> ( 9 BeX") 


+ 7 










fJrcc oc Z 2 B scaling 


8 B + 


x- 


" ( 8 BX-) 


+ 7 










CTree oc Z 2 B scaling 


ipx- 


-) + 


• a (al _ ) 






it 


estimated 


QM 3-body 


Kamimura 09 [75] 


{dX~ 


-) + 


- a -» (aX~) 


+ d 








QM 3-body 


Kamimura 09 [75] 


(tx- 


-) + 


- a — > (aX~) 


+ t 








QM 3-body 


Kamimura 09 |75| 



by colour, and dashed lines of the same color are the corresponding free states, e.g., the 
solid red line represents the ( 4 HeX~) abundance, while the dashed red represents that of 
free 4 He. Fig. [2] illustrates the case of a very long-lived X~. Focusing first on the bound 
states (solid colored lines) we see that 7 Be recombines first, followed by 7 Li, then 4 He, and 
finally protons. These results were to be expected, given that recombination occurs at 
Tfjec ~ B/\ln T)\ ~ 5/25, where B ~ Z 2 , with i] = n^/n^ the baryon-to-photon ratio. 

Comparing the solid lines, we see that free X~ particles dominate until 4 He recombines, 
leading to the first kink in the black curve, after which most X~ are in ( 4 HeX~) bound 
states q Subsequently the protons recombine, and this leads to the second kink in the 
black free X~ curve, as well as a small rise in the 4 He curve. This is because, immediately 
after the {pX~) states recombine, the important (pX~) + 4 He — > p + ( 4 HeX~) charge- 
exchange process converts (pX~) states into more ( 4 HeX~) bound states. The (dX~) state 
actually form earlier than (pX~) states, because the deuteron states are more tightly bound. 
However, deuterons are much rarer, and thus the (dX~) abundance is always quite small, 
and ultimately the (dX~)/(pX~) ratio is comparable to the ordinary D/H ratio. 

Turning attention now to the A = 7 states (blue and green curves), we see that recombi- 
nation into these nuclei occurs mostly after they are formed. The ( 7 BeX~) state has almost 
the same abundance as the free 7 Be state, whereas the abundance of ( 7 LiX~) is smaller 
than that of free 7 Li. At late times T ~ 10~ 5 MeV, the ( 7 BeX~) captures an electron and 
converts to ( 7 LiX _ ), which remains bound because here Q is smaller than the difference in 
binding energies. Looking at the 6 Li abundances, we see that the bound state has a much 
smaller abundance than the free nucleus. 

2 Note that it is important for this analysis that Yx- < ^4jJe' so a ^ ^~ P ar ^i c l es can fi n d 4 He 
partners. 
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Table 3: Summaries of the treatments of bound-state nuclear rates assumed in our work and 
in the recent literature (II) 
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simple scaling 
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Kamimura 09 [75] 
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Kamimura 09 [75] 
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simple scaling 
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Kamimura 09 [75] 


CtX~"| 4- 3 He 


-4 7 Bc + X - 


simple scaling 






QM 3-body 


Kamimura 09 [75] 


(aX~) 4- 4 Hp — 


V ( 8 BeX~) 4- 'v 










Pospelov 07 [76] 


(rvT _ "l 4- 6 T,i — 


, flORX - ! 4- -Y 










scaling from Caughlan 88 82 


(rvX~\ 4- 7 T,i — 


. n Rr"l -1- *v 










scaling from Angulo99 [53] 


(rv X~'\ -4- 7 "Rp — 


r ^ V . \ J T | 










scaling from Angulo99 [83] 


f 6 T,i X~ \ 4- 77 — ^ 


(T 4- 3 He 4- X~ 


simple scaling 






QM 3-body 


Kamimura 09 [75] 


( 6 LiX-) + n - 


±t + ct + X~ 










Caughlan 88 [82] 


( 6 LiX-) + d -4 


7 Li + p + X- 










scaling Malaney 89 [84] 


( 6 LiX-) + 


7 Bc + n + X- 










scaling Malaney 89 [5J] 


(pX-)4- 6 Li^ 4 


He + 3 Hc + X~ 






estimated 




Pospelov [36] 


( 7 ux-)+ P - 


> a + a + X~ 


simple scaling 






QM 3-body 


Kamimura 09 [75j 


( P X-) + 7 Li - 


-4 8 Bc + X- 








QM 3-body 


Kamimura 09 [75] 


{ 7 BeX-) +n - 


> 7 u + p + x- 










scaling from Cyburt04 [85] 


{ 7 BcX-) +p 


4 8 B + I- 








QM 3-bodyt 


Kamimura 09 [75] 


{pX-) + 7 Be 


4 8 B + r 








QM 3-body 


Kamimura 09 [75] 


( r BeX-) +d- 


> p + 2a + X~ 










Caughlan 88 [5J 


( 8 BeX-) + n 


-4 9 Be + X- 






estimated 


* 


Pospelov [36] 


( 8 BeX-) + d 


-4 10 B + X- 










scaled from Coc 12 [56] 


( 8 BeX-)+d^ 


6 Li4-a + X- 










scaled from Coc 12 [86] 


( 8 B1") -4 


8 Be + X- 










/3 lifetime Matt 64 [87] 



t This rate is mx -dependent. 

* The reaction ( 8 BcX~) + n is argued in |75| / to be non-resonant, which would reduce the 9 Be production 
from the levels given by this rate. 



Finally, we turn to the special case of ( 8 BeX~). Because this state has no analogue 
nuclide in standard BBN, there is no 8 Be + X~ recombination, and thus ( 8 BeX~) does 
not emerge when the temperature drops below its binding. Rather, production occurs via 
( 4 HeX~) + 4 He — > ( 8 BeX~) + 7, and thus we see that the abundances rises after that of 
( 4 HeX-). 

Our results are similar to those shown in [36J, apart from 6 Li, where we have purposely 
removed catalysis effects for this exercise only. One difference is that we use the charge- 
exchange rates of [75], which are larger than those used in [36], and thus give a smaller 
(pX~) abundance due to the more efficient {pX~) conversion. 

Figs. [3] and H] are for Tx = 10 6 and 10 3 s, respectively, and reveal few surprises. They are 
similar to the case of long-lived X until the time t ~ Tx, where t ~ l/T]^ eV , at which stage 
the X particles decay away. A corollary of this exercise is that if Tx < 10 3 s, recombination 
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10 1 10° 10 1 10 2 10" 3 10" 4 10" 5 
temperature T [MeV] 

Figure 2: The abundances of free nuclei (dashed lines) and nuclear bound states (solid 
lines) as functions of temperature. The black line corresponding to the abundance of free X~ 
particles, which is assumed to be 10~ 2 initially. In this case the X lifetime Tx is assumed to 
be infinite. 

cannot form bound states before the X particles decay, and hence bound-state chemical 
effects are negligible. 

5 Decays of the NLSP in Gravitino Dark Matter Sce- 
narios 

As in we base our discussion here on the CMSSM, as described in Appendix B. The mass 
of the gravitino is a free parameter in the CMSSM, and is the LSP and constitutes the dark 
matter if its mass, m 3 / 2 is chosen to be less than min(m x ,m f ). The abundances of the light 
elements provide some of the most important constraints on such a gravitino dark matter 
(GDM) scenario [ITl[l9l[2rJl[22l[26l[Ml[35l[37l[88HM ■ Tneir abundances also impose important 
constraints on neutralino LSP scenarios, since a gravitino NLSP could decay sufficiently 
slowly to affect them. Here, however, we restrict our attention to GDM scenarios with either 
a stau or neutralino NLSP, later focusing more closely on the stau NLSP case. 

In the process of calculating the lifetime of the NLSP, we calculate the partial widths 
of the dominant relevant decay channels of the NLSP and hence the various NLSP decay 
branching ratios. We also calculate the resulting electromagnetic (EM) and hadronic spectra, 
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10 1 10° 10 1 10 2 10" 3 10'' 
temperature T [MeV] 



10" ; 



Figure 3: .4s in Fig. 0, but with an assumed lifetime Tx = 10 6 s. 




10" 1 10" z 10" 
temperature T [MeV] 



Figure 4: .4s in Fig. [H but with an assumed lifetime tx = 10 3 s. 
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which impact the light-element abundances. The decay products that yield EM energy 
obviously include directly-produced photons, and also indirectly-produced photons (e.g., via 
the decays of neutral pions, 7T°), and charged leptons (electrons and muons) that may be 
produced via the secondary decays of gauge and Higgs bosons. Hadrons (nucleons and 
mesons such as the K^, K ± and 7r ± ) are produced through quark-antiquark pairs and via 
the secondary decays of gauge and Higgs bosons, as well as (in the case of mesons) via 
the decays of r leptons. We note that mesons decay before interacting with the hadronic 
background [13, 88, 89J . They are therefore irrelevant for BBN processes and to our analysis, 
except via their decays into photons and charged leptons. Thus the hadronic injections on 
which we focus our attention are those that produce nucleons, namely the decays via gauge 
and Higgs bosons and quark-antiquark pairs. 

In the case of the neutralino NLSP, we include the two-body decay channels x ~~ G Hi 
and x ~^ GV, where Hi = h,H,A and V = 7, Z, and also the dominant three-body decays 
X ->■ G7* ->■ Gqq and x -> GW + W~. In the case of % ~* GW + W~ we have included 
all the contributing tree- level amplitudes, as was done in [ID], thus treating correctly the 
longitudinal components of the W bosons. In general, the two-body channel x ~ * G^y 
dominates the x NLSP decays and yields the bulk of the injected EM energy. When the x 
is heavy enough to produce a real Z boson, the next most important channel is x ~ >* G Z, 
which is also the dominant channel for producing hadronic injections in this case. The Higgs 
boson channels are smaller by a few orders of magnitude, and those to heavy Higgs bosons 
(H,A), in particular, become kinematically accessible only for heavy x i n the large-m^ 
region. 

Turning to the three-body channels, the decay through the virtual photon to a qq pair 
can become comparable to the subdominant channel x ~ * G Z, injecting nucleons even in the 
kinematical region m x < m 3 /2+M z , where direct on-shell Z-boson production is not possible. 
In principle, one should also include qq pair production through the virtual Z-boson channel 
X — >■ G Z* G qq [18J and the corresponding interference term. However, this process is 
suppressed by a factor of M|/m(gg) 4 with respect to x — > G 7* — > G qq, and the interference 
term is suppressed by M^/m(qq) 2 . Numerically, these contributions are unimportant, and 
we drop them in our calculation. Finally, we note that the three-body decays to W + W~ 
pairs and a gravitino are usually at least five orders of magnitude smaller. Analytical results 
for the amplitudes for these gravitational decays of a neutralino NLSP have been presented 
in jlQ]. There, they were calculated for the inverse processes G — > x + X> but the decay 
amplitudes are the same, and the only adjustment needed is to interchange the neutralino 
and gravitino mass in the phase space. 

We presented in jlQ] our method of estimating the EM and the hadronic decays of the 
direct products of the x decays using the PYTHIA event generator [92]. We first generate a 
sufficient number of spectra for the secondary decays of the gauge and Higgs bosons and the 
quark pairs. Then, we perform fits to obtain the relation between the energy of the decaying 
particle and the quantity that characterizes the hadronic spectrum, namely dNh/dEh, the 
number of produced nucleons as a function of the nucleon energy. These spectra and the 
fraction of the energy of the decaying particle that is injected as EM energy are then used 
to calculate the light-element abundances. We use the same approach here. 
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An analogous procedure is followed for the f NLSP case. In [26], we assumed that the 
lighter stau was right-handed, so we ignored the stau mixing effects and the stau interactions 
with the W . However, in this analysis here we include the full effects of stau mixing. The 
decay rate for the dominant two-body decay channel, namely f — > Gt, was given in [90J. 
However, this decay channel does not yield any nucleons. Therefore, one must calculate some 
three-body decays of the f to obtain any protons or neutrons. The most relevant channels are 
f ^GrMGZr, f ^ Zt* ^ G Zr, f rx°* ->• G Z r and f -> G Z In addition, so 
as to include the full effects of stau mixing, we included the processes f — > G r* — > G W~ u T , 
f -> W~ z7* -> GW~v r , f -> v T xT ->■ GW~u r and f ->■ GW~ v T . 

Analytical results for three-body stau decays can be found in Appendix C. We then use 
PYTHIA to obtain the hadronic spectra and the EM energy injected by the secondary W, 
Z-boson and r-lepton decays. As in the case of the \ NLSP, this information is then used 
for the BBN calculation. 

6 NLSP Lifetimes in the CMSSM with Gravitino LSP 

As discussed above, we study the constraints from the cosmological light element abun- 
dances in the context of the CMSSM. The recent discovery of a new boson with mass ~ 125 
to 126 GeV with properties that resemble those of the Standard Model Higgs boson [93] 
motivates us to concentrate on regions of the CMSSM parameter space where the lightest 
neutral Higgs boson has a mass close to this range, taking into account the theoretical un- 
certainty in the calculation of its mass for any fixed values of the CMSSM parameters [94|. 
As discussed in [95], this mass range favours large values of A and tan/3: see also [96|l9"T]. 
On the other hand, the constraint from B s — > [i + [i~ [98] disfavours very large tan /3 [99, 100J. 
Accordingly, in this paper we discuss one example of a (mi/ 2 ,m ) plane with tan/3 = 10 and 
two examples with tan/3 = 40 [93]. In many models of supersymmetry breaking, the soft 
trilinear supersymmetry-breaking parameter Aq oc mo- For tan /3 = 10 we consider the single 
value A = 2.5 m , and in the tan /? = 40 case we consider the two options A = 2 m , 2.5 mo- 
In the absence of clear indications on the gravitino mass 1/13/2, in each case we consider two 
options: fixed 7713/2 = 100 GeV and 7773/2 = 0.1m . We also consider in less detail two 
examples of (mi/ 2 ,Ao) planes with fixed tan/3 = 40 and m 3 / 2 = 100 GeV, namely with 
m = 1000, 3000 GeV. 

An important ingredient in understanding the morphology of our results in the (mi/2, mo) 
planes is provided by the NLSP lifetime t NL sp- Fig. |5] displays contours of tnlsp i n the 
first cases studied above, namely the (m^jWo) planes for tan/3 = 10, A = 2.5 m and 
?7i3/2 = 100 GeV (left) and 777,3/2 = 0.1m (right). In the upper part of the left panel, 
where the lightest neutralino is the NLSP, we see that the lifetime contours are essentially 
vertical, since they depend mainly on the relationship of m x (and hence mi/2) to 7773/2- 
These contours appear only when the gravitino is the LSP, i.e., for m x > 7773/2, and there is 
a vertical band at small 777^/2 where this condition is not satisfied. Also, we note that the 
lighter stop squark is either the NLSP or tachyonic in the grey shaded triangular regions in 
the small- 777i/2, large- mo corners of these panels. In the lower part of the left panel of Fig. 
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where the lighter stau is the NLSP, the contours of constant NLSP lifetime curve, track the 
relationship between and m fl . The lifetime contours in the right panel of Fig. [5j for 

^3/2 = 0.1 mo, are everywhere sloping up from left to right. 




m V2 [TeV] m m [TeV] 

Figure 5: The NLSP lifetime tnlsp in the (mi/2,mo) plane for A = 2.5m , tan/3 = 10 
and 1713/2 = 100 GeV (left) and m 3 / 2 = 0.1 mo (right). 

Fig. [6] displays the corresponding contours of tnlsp for the cases tan/3 = 40, A = 2 mo 
and m 3 / 2 = 100 GeV (left) and m 3 / 2 = 0.1 m (right). These exhibit similar features to the 
previous case, except that the stau NLSP region is now larg result of the larger value 

of tan/3, and now we see a difference in the behaviours of the lifetime contours in the stau 
and neutralino NLSP regions. The vertical band at small mu2 where the gravitino is not 
the LSP is now fully visible. In this case, there is a triangular region at small and large 
m where the gravitino is no longer the LSP. 

Finally, Fig. [7] displays the corresponding contours of t NL sp for the cases tan /3 = 40, 
Aq = 2.5 m and m 3 / 2 = 100 GeV (left) and m 3 / 2 = 0.1 m (right). This is qualitatively 
similar to Fig. [6], though we note that the stau NLSP region has expanded again, this time 
as a result of the larger value of Aq. Note also that the triangular region where the stop is 
light (or tachyonic) has reappeared at the larger value of A . 
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Figure 6: The NLSP lifetime in the (mi/2,mo) plane for A = 2m , tan/3 = 40 and 
m 3/2 = 100 GeV (left) andm 3 /2 = 0.1 m (right). 

7 Light-Element Constraints in the CMSSM with a 
Metastable Stau NLSP 

We display in Fig. [8] the light-element abundances we calculate in the (m 1 /2,m ) plane for 
the first example introduced above, namely Aq = 2.5 mo, tan = 10 and 1713/2 = 100 GeV. 
In this and subsequent figures, the stau is the NLSP in a wedge of each plane at low itlq 
and large m\/2- (The outline of this wedge can be seen in each of the preceding lifetime 
plots by connecting the points where the lifetime vs. mi/2 changes from a curve to a straight 
line.) In most of the planes at larger mo the lightest neutralino is the NLSP. However, 
when Aq = 2.5mo, there are also wedges at large m and small mi/2, shaded grey, in which 
the NLSP is the lighter stop squark. Indeed for very low m^, the stop mass squared is 
negative and hence for parameter choices inside this grey wedge the sparticle spectrum is 
not physical. There is only a very narrow strip along the wedge where the stop is actually 
the NLSP. We do not consider this case in the present work (see however |101j ). discussing 
only the neutralino and stau NLSP cases. 

7.1 Summary of Light-Element Abundance Constraints 

As in subsequent similar figures, the upper left panel of Fig. displays the D/H ratio, 
the upper middle panel displays the 3 He/D ratio, the upper right panel displays the 4 He 
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% [TeV] m m [TeV] 

Figure 7: TTie NLSP lifetime in the (mi/ 2 ,mo) p/ane /or A = 2.5 mo, tan/3 = 40 and 
m 3/2 = 100 GeF (left) andm 3 /2 = 0.1 m (right). 

abundance, the lower left panel the 6 Li/ 7 Li ratio, the lower middle panel the 7 Li/H ratio, 
and the lower right panel the 9 Be/H ratio. As a general rule, we consider the regions left 
unshaded to be compatible with observation, whereas the yellow regions are problematic, 
and the red and magenta regions are progressively more strongly excluded. Solid shadings 
are used for regions with excess abundances, and hashed shadings for regions with low 
abundances. The criteria adopted for the light-element abundances are similar to those used 
in our previous work, and are summarized in Table H||E 



D/H 

We assume the mean value given in |102j 

J = (3.0 ±0.7) x 10^ 5 , (14) 
/ p 

corresponding to the deuterium abundance measured in 10 quasar absorption systems [103J, 
and the quoted uncertainty is given by the sample variance in the data. This is considerably 
larger than the error in the mean, which is only 0.2. Therefore, we consider any value outside 
the range (2.3 — 3.7) x 10~ 5 as problematic, as indicated in Table HI which also includes ranges 
that we consider to be (strongly) excluded. 

3 The values corresponding to 'strong exclusion' are somewhat arbitrary, but serve to indicate how rapidly 
the abundances are varying in relevant regions of parameter space. 
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tan|3=10 m 3/2 = 100 Ge V A = 2.5 m 
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Figure 8: Light-element abundances in the (mi/ 2 ,m ) plane for A = 2.5m , tan/3 = 10 
and m 3 /2 = 100 GeV. 

3 He/D 

Whilst it is difficult to use 3 He to constrain BBN, it is possible to use the ratio 3 He/D [104J. 
Although 3 He may be created or destroyed in stars, D is always destroyed in the pre-main 
sequence of stellar evolution and, as a result, the ratio 3 He/D is a monotonically increasing 
function of time. Thus one can use the solar ratio of about 1 [105J to constrain the BBN 
ratio. Because 3 He can be produced and/or D can be destroyed, we do not assume a lower 
bound to the ratio. 



4 He 

Although the determination of the 4 He abundance in extragalactic HII regions is dominated 
by systematic uncertainties |106j . using the Markov Chain- Monte Carlo methods described 
in |107j and data compiled in [108], one finds |109] 

Y v = 0.2534 ± 0.0083 (15) 
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Comparison with 


D/H 


3 He/D 


4 He 


b Li/ 7 Li 


7 Li/H 


9 Be/H 


observation 


(xlO- 5 ) 








xlO" 10 


xlO" 13 


Strongly excluded 


< 0.5 


— 


< 0.22 


— 


< 0.1 


— 


Excluded 


< 1.0 




< 0.23 




< 0.2 




Problematic 


< 2.3 




< 0.24 




< 0.5 




Acceptable 


[2.3, 3.7] 


[0.3, 1.0] 


[0.24, 0.27] 


< 0.05 


[0.5, 2.75] 


< 0.3 


Problematic 


> 3.7 


> 1.0 


> 0.27 


> 0.05 


> 2.75 


> 0.3 


Excluded 


> 5.0 


> 3.0 


> 0.28 


> 0.1 


> 10 


> 1.0 


Strongly excluded 


> 10 


> 5.0 


> 0.29 


> 0.2 


> 30 


> 3.0 



Table 4: The ranges of light-element abundances whose comparisons with observation we 
consider in this work to be acceptable, problematic and (strongly) excluded, as shown in the 
unshaded/yellow/red/magenta regions in the Figures. 



based on a regression of Y vs. O/H and 

(Y) = 0.2574 ± 0.0036 (16) 

based on a weighted mean. As we will see, once the standard BBN value of Y p is affected by 
NLSP decays, it varies very rapidly and it suffices to consider values outside the range [0.24, 
0.27] to be problematic. 

6 Li/ 7 Li 

Some observations of Li absorption lines in halo dwarf stars have claimed evidence for a 
relatively large ratio of 6 Li/ 7 Li ~ 0.05 jl 10] over a broad range of metallicities, though it 
remains possible that these observations are also dominated by systematic uncertainties [11 lj . 
There are a few reliable observations of stars with a similar ratio of 6 Li/ 7 Li in a very narrow 
range of metallicity jM] consistent with galactic cosmic-ray nucleosynthesis |112j . However, 
no observations indicate a ratio greater than 0.05 which we set as our lower boundary of the 
problematic range. 



7 Li/H 

The cosmological 7 Li problem [52] is now well established. There are many observations 
of 7 Li in halo dwarf stars [63] that indicate a far lower 7 Li/H abundance than predicted in 
standard BBN. We adopt the range found in the plateau of Lithium versus metallicity [113], 
namely 

^ = (1-23^) x 10- 10 , (17) 

n / halo* 

although the lithium abundance observed in globular cluster stars may be a factor ~ 2 
higher [114] . Although the preferred range in ( 1T7|) is rather narrow, we deem that any 
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reduction from the BBN value of 7 Li/H = (5.07±g;S) x 10~ 10 [52] to < 2.75 x KT 10 
represents a significant improvement in the 7 Li problem, and we take this to be the lower 
bound of our problematic region. NLSP decays can also destroy too much 7 Li and we will 
consider any value below 0.5 x 10~ 10 as similarly problematic. 

9 Be/H 

Finally, 9 Be is also observed in halo dwarf stars, and is found to vary strongly with metallicity, 
as seen in a recent set of observations [115j . These observations extend down to [O/H] of 
about -2.5 with a 9 Be/H abundance of 3 x 10 -14 . Though there is a single observation [116] 
of 9 Be with an abundance about 3 times lower, conservatively we will consider problematic 
any 'primordial' abundance in excess of the highest value seen at the lowest metallicity. 

In the cases of 3 He/D, 6 Li/ 7 Li and 9 Be/H, there are no observational lower limits, so we 
do not quote ranges of abundances that we consider too low. Within the unshaded regions, 
we also display extra contours for the 4 He abundance = 0.25 and 0.26 (dashed and solid lines, 
respectively), for 7 Li/H = 0.91,1.91 x 10" 10 (dashed and solid) and 9 Be/H = 1,2 x 10~ 14 
(dashed and solid). In the case of 7 Li/H, as already discussed, it is well known that standard 
BBN gives a ratio significantly higher than that indicated by observations. Therefore, large 
parts of the regions coloured yellow in the 7 Li panels yield an abundance that is no worse than 
in standard cosmology, and may even be in somewhat better agreement with observation. 
Depending how seriously one takes the cosmological 7 Li problem, the favoured (unshaded) 
regions in subsequent plots could be expanded. In general, we see discontinuities in the 
colouring along a rising diagonal line: above it, the lightest neutralino is the NLSP, and 
below it the lighter stau is the NLSP, which is the case of main interest here. 

We recall from previous analyses that hadronic processes are mostly relevant for lifetimes 
^ 10 4 s, whereas electromagnetic processes are generally dominant for longer lifetimes, i.e., at 
smaller for any fixed value of m . We also note that, for any fixed mi/2, the abundance 
of metastable relic particles (before decay) is generally largest at large m . We therefore 
expect hadronic processes to be most important when both m\/2 and mo are large. Indeed, 
in the upper right parts of the 4 He panels in Fig. [8] and later figures we see triangular regions 
where the 4 He abundance is enhanced unacceptably by hadroproduction. This is generally 
accompanied by hadronic depletion of the 3 He/D ratio and enhancements in D/H, 6 Li/ 7 Li, 
7 Li/H and 9 Be/H. On the other hand, staying at large m , the dominant electromagnetic 
processes at smaller m\/2 include photodestruction of 4 He and 7 Li, accompanied by photo- 
production of 3 He/D and D/H. 

7.2 Application to the CMSSM with a Metastable Stau NLSP 

We see in the upper panels of Fig. [8] that for A = 2.5 mo, tan/3 = 10 and m.3/2 = 100 GeV 
the D/H ratio is acceptable in arcs with m\/2 ~ 2 TeV and > 4 TeV, whereas the 3 He/D ratio 
is generally acceptable for my 2 > 1-6 TeV and the 4 He abundance is acceptable throughout 

4 This corresponds to the BBN value at a baryon-to-photon ratio of -q — 6.16 x 10 -10 [50]. A similar value 
of (5.24 ± 0.5) x 10 -10 was found in another recent analysis |86) . 
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the stau NLSP wedge of the (7711/2, tuq) plane (This demarkation is displayed in the summary 
plot below). In the lower panels of Fig. [8] we see that the 6 Li/ 7 Li ratio is unacceptable for 
TO1/2 < 3 TeV, that the 7 Li/H ratio is acceptable in an arc with rai/2 > 2.5 TeV, as well as in 
the neutralino NLSP region with mi/2 > 4 TeV at large mo- The 9 Be/H ratio favours either 
mi/a > 3 TeV or a triangular region with a neutralino NLSP with myz < 1 TeV. The overall 
conclusion is that all the light-element abundances are acceptable in a narrow arc starting 
on the stau/neutralino NLSP boundary where (m 1 / 2 ,m ) ~ (4.0, 1.1) TeV and extending to 
larger mi/2 at lower uiq 0. Note that there is excessive photo-destruction of both 4 He and 
7 Li in the low mass neutralino NLSP region. While there is a narrow strip along 777.1/2 ~ 1.4 
TeV where 7 Li is just right, this strip is excluded by many of the other light elements. This 
behaviour is also seen in the subsequent parameter choices. 

Fig. |9] displays a similar analysis for the same CMSSM parameters A = 2.5 m , tan/3 = 
10, but with m 3 / 2 = 0.1 m . In this case, we see that the light-element abundances are all 
acceptable in a narrow arc through the stau NLSP region between (7771/2, m o) ~ (2-3, 0.4) TeV 
and ~ (4.8, 1.5) TeV, that is defined essentially by the D/H and 7 Li/H constraints 0. 

Fig. [TU1 summarizes our results for the CMSSM (mi/2, itiq) planes for Aq = 2.5mo, tan /3 = 
10, with 7773/2 = 100 GeV (left) and 772,3/2 = 0.1 mo (right). Here and in subsequent similar 
figures, the magenta regions are strongly excluded by one or more constraints, the red regions 
are excluded by one or more constraints, and the yellow regions are problematic for at least 
one constraint. We see explicitly the unshaded narrow arcs where all the constraints are 
satisfied. These are the regions where the cosmological Li problem is solved in the presence 
of metastable stau NLSPs: they are all below the grey line that marks the boundary between 
the neutralino and stau NLSP wedges. 

Also shown is a green line, above which the gravitinos arising from NLSP decays have 
a density greater than the range allowed by WMAP and other observations jSU]- This line 
corresponds to the gravitino relic abundance determined from NLSP decays 



where I^nlsp is the thermal relic density of NLSPs left over after annihilations. We note 
that there may be other sources of gravitinos such as reheating after inflation which would 
further strengthen this bound. This constraint excludes almost completely the neutralino 
NLSP regions in Fig.[T0]and the subsequent analogous figures, but does not impact the white 
regions compatible with all the light-element constraints. Also shown in this and subsequent 
summary figures are some contours of calculated values of = 124 GeV (dash-dotted), 
125 GeV (solid), 126 GeV (dotted) and 127 GeV (dashed). The present experimental and 
theoretical uncertainties are such that no calculated value of Mh G [124, 127] GeV can 
currently be excluded, and an even larger range may be permitted at large 7771/2, where the 
FeynHiggs code [91] warns of theoretical uncertainties considerably exceeding 1.5 GeV. 

5 We note in passing that there is no region of the neutralino NLSP wedge where all the light-element 
abundances are acceptable. 

6 Again, there is no acceptable region where the neutralino is the NLSP. 




(18) 
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Figure 9: Light-element abundances in the (mi/ 2 ,mo) plane for A = 2.5 m , tan/3 = 10 
and = 0.1 m . 

Looking back at the contours of constant tnlsp i n described arcs in Fig. [U we see that 
in the stau NLSP segment of the (mi/2,wio) plane they parallel the contours in the corre- 
sponding regions of the (mi/2, mo) planes in Figs. El and [9] for the different light-element 
abundances. This confirms the important influence of 7nlsp- Comparing with the summary 
of this case displayed in the left panel of Fig. [TDl we see that in this case the optimal lifetime 
for solving the cosmological Lithium problem is tnlsp ~ f ew x 10 2 s. In the case when 
m 3/2 = 0.1 mo (right panel of Fig. Ej), we see that in the stau NLSP segment the contours 
of constant Tnlsp parallel those of constant 6 Li/ 7 Li and 9 Be/H ratios, though the shapes of 
the D/H and 7 Li/H contours are quite different. Looking at the right panel of Fig. [TO], we 
see that in this case the optimal lifetime for solving the cosmological Lithium problem is 
tnlsp ~ few x 10 3 s. 

Though we do not show the results here, we have studied other choices for the gravitino 
mass for the values of tan/3 = 10 and Aq = 2.5 mo. For example, for a larger fixed gravitino 
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Figure 10: Summary of the light- element- abundance constraints in the (mi/2,mo) plane for 
A = 2.5 mo, tan/3 = 10 and m 3 /2 = 100 GeV (left) and 777.3/2 = 0.1 mo (right). 

mass of 500 GeV, one must consider larger m.1/2 > 1 TeV to ensure a gravitino LSP. For 
a given gaugino mass, the NLSP lifetime is longer. As a result, the acceptable arc of D/H 
moves to larger 7711/2- More importantly, the 6 Li constraint would now exclude all values of 
mi/2 between 1 and 5 TeV. The 9 Be constraint similarly would exclude the entire stau NLSP 
region displayed. Had we chosen instead 777.3/2 = mo, a gravitino LSP would be present only 
in the lower right half of the plane. Once again lifetimes would typically be longer, affecting 
the light element abundances. In this case, only a small corner of the parameter space at 
very large mi/2 and very small mo would survive all constraints. 

We now describe an analogous analysis for the CMSSM (mu2, m ) planes for A = 2m , 
tan/3 = 40. Fig. [TT] displays our results for the option m 3 / 2 = 100 GeV. In this case, the 
D/H constraint would allow most of the lower half of the parameter plane. This regions 
would be allowed by both the 3 He/D (except for a small area with low m\/% and mo) and 
4 He constraints, but much of it is excluded by the 6 Li/ 7 Li ratio, and more strongly excluded 
by the 9 Be/H ratio. Improvement in the 7 Li/H ratio only occurs around an arc starting at 
(m 1 / 2 ,m ) = (3.2,2) TeV. This arc is for the most part allowed by the other constraints. 

Fig. [T21 displays the results of a similar analysis for 777.3/2 = 0.1 m , but with the same 
values of the CMSSM parameters. Once again, the neutralino NLSP region is excluded by 
the D/H ratio, which is also problematic for a large area with mo > 1 TeV. The 3 He/D 
and 4 He constraints are qualitatively similar to the previous case. However, the effect of 
the 6 Li/ 7 Li constraint is somewhat different: it excludes a bulbous region of the stau NLSP 
segment extending almost to m 1( / 2 ~ 5 TeV as does the 9 Be constraint. In this case, the 
arc allowed by the Li/H ratio is wider and has shifted to larger masses. As a result, the 
only region that has a chance of being compatible with all the light-element constraints has 
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Figure 11: Light-element abundances in the (mi/ 2 ,mo) p/ane for A = 2m 0; tan/3 = 40 
and m 3 /2 = 100 GeV. 

mi/2 ~ 5 TeV and mo ~ 1 TeV. 

The results in the (mi/2, mo) planes for v4o = 2 mo and tan/3 = 40 are summarized 
in Fig. [13j In the case of = 100 GeV (left panel) we see an allowed arc across the 
stau NLSP region extending from (m 1 /2,m ) ~ (3.2,2) TeV to ~ (5,1.5) TeV. In the case 
of 777.3/2 = 0.1 mo (right panel), there is only a very small region of marginal consistency 
close to (mi/2,mo) ~ (5, 1) TeV. In the stau NLSP region the contours of the light-element 
constraints displayed in Figs. [TT]and[T2] (except for the D/H ratio in the latter case) again 
parallel the contours of constant r NL sp in Fig. [HI and we see in Fig. [13] that the preferred 
ranges of tnlsp are somewhat below and above 10 3 s, respectively. 

For the same choice of tan/3 = 40 and Aq = 2.0mo, had we taken m 3 / 2 = 500 GeV we 
would have found that, due to the increased lifetimes, virtually the entire parameter plane 
with a gravitino LSP would be strongly excluded by the 6 Li/ Li ratio. The 9 Be constraint 
also would strongly exclude the stau NLSP region shown. For m 3 / 2 = mo, we would once 
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Figure 12: Light-element abundances in the (mi/ 2 ,mo) p/ane /or v4 = 2m , tan/3 = 40 
and 777.3/2 = 0.1 m . 

again be forced into a tiny area in the lower right corner of the (mi/2, mo) plane. 

Turning now to the case Aq = 2.5 mo, tan /3 = 40 and 7773/2 = 100 GeV shown in Fig. [H} 
we note in particular that there is virtually no improvement over standard BBN in the 7 Li/H 
abundance throughout almost all the stau NLSP region. Only a small region close to the 
stau-neutralino NLSP boundary extending to higher masses from ~ 3 TeV is consistent 
with this constraint. The 6 Li/ 7 Li and 9 Be constraints once again dominate in the stau NLSP 
region. 

Even this small region of consistency is eradicated in the case A = 2.5 777 , tan/3 = 40 
and 7773/2 = O.lmo shown in Fig. [T5j Now, the improvement in 7 Li/H is limited to a small 
region at very large masses (mi/ 2 , 777 ) ~ (5,3) TeV, and this corner of parameter space is 
excluded by both the 6 Li/ 7 Li and 9 Be/H ratios. 

These results are summarized in Fig. [16j We see in the left panel for 7773/2 = 100 GeV 
that only a very small region with (777.1/2,777.0) ~ (3.5,2.8) TeV is compatible with all the 
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Figure 13: Summary of the light- element- abundance constraints in the (mi/2,mo) plane for 
A Q = 2m , tan/3 = 40 and m 3 / 2 = 100 GeV (left) and m 3 / 2 = 0.1 m (right). 

constraints, whereas we see no allowed region in the right panel for m 3 / 2 = 0.1 mo- Comparing 
the lifetime contours in Fig. [7]with Figs. HH [15]and[l6l we again see that the preferred ranges 
of tnlsp are ~ 10 3 s. Choosing the gravitino masses m 3 / 2 = 500 GeV or = m would leave 
us with results very similar to those described for A = 2m . 

The results shown above have been for slices through the CMSSM parameter space 
corresponding to (mi/ 2 ,m ) planes for fixed tan/3 and A . We have also explored how the 
results for tan /3 = 40 vary as functions of A for a couple of values of m = 1000, 3000 GeV, 
with the results summarized in Fig. (T7] The left panel is for m = 1000 GeV, which is 
typical of the range of m in the unshaded regions in the cases studied above. We see that 
a large region with mu2 > 4 TeV and Aq < 2 TeV is unshaded and hence 7 Li-compatible. 
On the other hand, we see no unshaded region in the right panel for itlq = 3000 GeV, which 
is less typical of the values of m found in the unshaded regions of previous summary plots. 
Therefore, we expect that the features found earlier are quite generic. 

Also shown in Figs. [10], [13], [16] and [T7] are some representative contours of the lightest 
MSSM Higgs boson M h , as calculated using the FeynHiggs code [94J. This code is generally 
thought to have an uncertainty ~ 1.5 GeV for generic sets of CMSSM parameters, but warns 
of larger uncertainties at the large values of mi/ 2 of interest here0. Accordingly, we consider 
calculated values of M h 6 [124, 127] GeV to be compatible with the observed range of 125 
to 126 GeV [93J, and an even larger range of calculated values of may be acceptable at 
large mi/ 2 - In the cases displayed in Fig. [10] we see that the ends of the BBN-compatible 

7 This may be linked with the irregular behaviours of some calculated contours of Mh in Figs. I131IT61 and 

El 
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Figure 14: Light-element abundances in the (mi/ 2 ,m ) plane for A = 2.5m , tan/3 = 40 
and m 3 /2 = 100 GeV. 

arcs with higher mo have Mh ~ 124 GeV, i.e., within the acceptable range, and hence may 
be preferred. In Figs. [13] we see that the preferred arc for 7713/2 = 100 GeV corresponds 
to Mh ~ 124 to 126 GeV, all within the range suggested by the LHC, whereas in the case 
m 3/2 = 0.1 Too the BBN-compatible region has Mh ~ 124 GeV. In Fig. [16] we see that 
the small BBN-compatible region for 777,3/2 = 100 GeV corresponds to a nominal value of 
Mh ~ 127 GeV, at the upper end of the LHC-compatible range. Finally, in Fig. [17] we see 
that the unshaded region in the left panel corresponds generally to M^ ~ 124 GeV, which is 
compatible within theoretical uncertainties with the LHC discovery. 

7.3 Bound-State Effects and Uncertainties 

We conclude this Section with a brief discussion of the importance of bound-state effects 
and their uncertainties. 
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Figure 15: Light-element abundances in the (mi/ 2 ,m ) plane for A = 2.5m , tan/3 = 40 

and 777.3/2 = 0.1 777 . 

Fig. [TBI shows how our results in the (7771/2,7770) plane for Aq = 2 mo, tan/3 = 40 and 
7773/2 = 100 GeV would change if all bound-state effects were to be switched off, but with 
decay effects retained. Comparing with Fig. [TTJ we see that the D/H, 3 He/D and 7 Li/H 
ratios are unaffected, as is the 4 He abundance. However, there are major changes in the 
6 Li/ 7 Li ratio and 9 Be abundance. In particular, arcs at small 777 1 / 2 and 777 that would have 
been permitted (modulo the lack of improvement in the 7 Li/H abundance) in the absence of 
bound-state effects are robustly excluded by both the 6 Li/ 7 Li and 9 Be/H ratios once bound- 
state effects are included. On the other hand, the allowed arc at larger 777 1 / 2 and mo is quite 
unaffected by bound-state effects, as seen by comparing the left panel of Fig. [19] with the 
left panel of Fig. [131 Even more dramatically, we see that the stau NLSP region excluded 
by 9 Be was entirely due to bound state effects. 

As discussed in Section 2, bound-state 9 Be production hinges on two principal uncertain- 
ties in our bound-state analysis. One of these uncertainties is the ( 8 BeX) binding energy, 
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Figure 16: Summary of the light- element- abundance constraints in the (m^^o) plane for 
A = 2.5 mo, tan/3 = 40 and m 3 / 2 = 100 GeV (left) and m.3/2 = 0.1 mo (right). 
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Figure 17: Summary of the light-element-abundance constraints in the (m.1/2, ^o) plane for 
tan/3 = 40 and m 3/2 = 100 GeV with m = 1000 GeV (left) and m = 3000 GeV (right). 

must be high enough to allow for ( 4 HeX~)+ 4 He — > ( 8 BeX~)+7 to be exothermic, i.e., Q > 
(cf. eq. [p. Our analysis has assumed as default the B 8 = 1.1679 value of ref . [75], which 
implies the formation reaction is strongly exothermic, and thus the reverse photodissociation 
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Figure 18: Light- element abundances in the (mi/2,mo) plane for A Q = 2m , tan (3 = 40 
and m 3 / 2 = 100 GeV, with decay effects retained but all bound-state effects switched off. 

of ( 8 BeX~) is strongly suppressed. We have also considered both greater and smaller values 
of the binding energy. Using the larger value B$ = 1.408 MeV [76] makes the reaction even 
more exothermic; this gives results that are almost identical to our default analysis, since the 
bound-state formation rate remains very similar. The central panel of Fig. [19] summarizes 
the overall effect on the allowed region of (^1/2,^0) plane for A Q = 2m , tan/3 = 40 and 
m 3/2 = 100 GeV, which is almost indistinguishable from the default result shown in the left 
panel of Fig. [I3J 

Our own three-body estimate of the ( 8 BeX~) binding energy in Section 2 gives B 8 = 
492 ± 50 keV, a value that exceeds the effective 'no-go' limit in ((T])only by Q = B$ — B™ n = 
53 keV. In this situation, ( 8 BeX~) production is weakly exothermic but remains highly 
vulnerable to photodissociation back to ( 4 HeX~) + 4 He. This reverse reaction suppresses 
( 8 BeX~) formation until the temperature drops to ~ Q/\ I1177I ~ 2 keV. But at this late 
time, no free neutrons are available, and the result is effectively that no bound-state 9 Be 
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Figure 19: Summary of the light-element-abundance constraints in the (mi/2,mo) plane 
for Aq = 2.5 m , tan/3 = 40 and m 3 / 2 = 100 GeV, if all bound-state effects were switched 
off (left), with greater 8 Be-X binding energy than our default choice (centre) and with 9 Be 
bound-state production suppressed (right). 

production occurs. Results for this case appear in Fig. [20] We see that in the lower-right 
region, where bound-state effects are important, the 9 Be production is now missing. Indeed, 
we have checked that the results are unchanged if bound-state 9 Be production is switched off 
entirely, as would be the case if the if ( 8 BeX) binding energy drops below the limit in eq. (JTJ). 
Note also that even in the absence of bound-state production, 9 Be contours do remain in the 
upper-left region in Fig. [201 In this regine the thermalized 4 He fragments deuterium, tritium 
and 3 He are overproduced, and these can still make 9 Be via reactions with background 7 Be. 

The right panel of Fig. [19] summarizes the overall effect of suppressed 9 Be bound-state 
production on the allowed region of the (mi/2,m ) plane for A = 2m , tan/3 = 40 and 
"^3/2 = 100 GeV. Perhaps surprisingly, there is little visible effect, and specifically none 
on the allowed unshaded arc at large m^j and itiq. This is because the 9 Be and 6 Li/ 7 Li 
constraint contours very closely shadow each other in the lower part of the (m.1/2, mo) plane. 

A second uncertainty in 9 Be production comes from the requirement that the ( 8 BeX~) + 
n — > ( 9 BeX~) + 7 reaction is on resonance with the first excited state of ( 9 Be*X~). As 
noted in section 2, this requires that the ( 8 BeX _ ) and ( 9 Be*X~) bindings conspire in such 
a way that the entrance channel is on resonance, as would be the case for the larger binding 
energy B% = 1.408 MeV. Because the ( 9 Be*X~) binding is quite uncertain, this possibility 
remains viable. However, if the excited state level turns out to fall far (> 100 keV) from 
the ( 8 BeX~) + n entrance, then the reaction will be non-resonant and suppressed. And here 
again, the bound-state 9 Be production would become unimportant, similar to the results in 
Fig. [20] and in the right panel of Fig. [191 

We conclude that, whereas the overall bound-state effects are very important, the prin- 
cipal uncertainties associated with the ( 8 BeX) binding energy have little effect on our final 
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Figure 20: Light- element abundances in the (mi/2,mo) plane for A Q = 2m , tan /3 = 40 
and 777.3/2 = 100 GeV, with bound-state 9 Be production suppressed. 

results. This comes about because the regions excluded by bound-state 9 Be production over- 
lap almost completely with those also excluded by bound-state 6 Li production. In particular, 
though our analysis incorporated the resonant ( 9 Be X~) reaction rate postulated in [76J, our 
final results are not very sensitive to this assumption. 

8 Summary 

We have presented in this paper a new treatment of the possible effects of bound states of 
metastable charged particles in the light-element abundances yielded by Big-Bang Nucle- 
osynthesis (BBN), including in our analysis calculations of the abundances of D, 3 He, 4 He, 
6 Li, 7 Li and 9 Be. We have applied our code to the case of metastable f\ NLSPs in the 
framework of the CMSSM with a gravitino LSP. Motivated by the discovery of (apparently) 
a Higgs boson weighing ~ 125 to 126 GeV, we have concentrated on regions of the CMSSM 
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in which this may be interpreted as the lightest neutral Higgs boson, specifically (mi/2,mo) 
planes with A > 2m . 

We find interesting examples in which the light-element abundances are as consistent 
with observations as are calculations of standard homogeneous BBN with no metastable 
relic particles. Indeed, we find generic strips of the CMSSM parameter space in which the 
cosmological 7 Li problem may be solved without altering unacceptably the abundances of 
the other light elements. Examples are given for tan/3 = 10 and 40 with Higgs masses 
compatible with the LHC discovery. 

Characteristics of these models include relatively large values of the soft supersymmetry- 
breaking parameters m 1 / 2 and m , with heavy supersymmetric particles that could not be 
detected directly at the LHC. Another characteristic of these models is that the f\ lifetime is 
(9(10 3 ) s. Avenues for future research include a more complete examination of the CMSSM 
parameter space, possible extensions to more general supersymmetric models, as well as to 
non-supersymmetric scenarios. In addition, we reiterate the need for careful study of the 
nuclear physics behind bound-state 8 Be and 9 Be production. 

The observational situation with the 6 Li and 7 Li abundances is still evolving, and the fat 
lady has not yet sung the final aria in the cosmological lithium saga. It is perhaps still possible 
that the current discrepancy with standard homogeneous BBN will eventually dissipate. 
However, we have shown that, should it survive, it could have a plausible supersymmetric 
solution. 

Acknowledgments 

B.D.F. is pleased to acknowledge illuminating discussions with Chris Hirata regarding bound- 
state recombination, and we thank Sven Heinemeyer for discussions about FeynHiggs calcu- 
lations. The work of R.H.C. was supported by the U.S. National Science Foundation Grants 
PHY-02-016783 and PHY-08-22648 (JINA). The work of J.E. and F.L. was supported in 
part by the London Centre for Terauniverse Studies (LCTS), using funding from the Eu- 
ropean Research Council via the Advanced Investigator Grant 267352: this also supported 
visits by K.A.O. and V.C.S. to the CERN TH Division, which they thank for its hospitality. 
The work of B.D.F. was partially supported by the U.S. National Science Foundation Grant 
PHY-1214082. The work of F.L. and K.A.O. was supported in part by DOE grant DE- 
FG02-94ER-40823 at the University of Minnesota, and the work of F.L. was also supported 
in part by the Doctoral Dissertation Fellowship at the University of Minnesota. The work of 
V.C.S. was supported by Marie Curie International Reintegration grant SUSYDM-PHEN, 
MIRG-CT-2007-203189. 

A Recombination of X Bound States 

A comparison of (AX~) recombination with that of ordinary hydrogen recombination reveals 
important similarities but also crucial differences. The recombination of ordinary cosmic 
hydrogen and helium does not proceed primarily through pe — > (pe)i transitions directly to 
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the n = 1 ground state. This is because such recombinations emit Lyman limit photons 
with energy = 5( pe ) = 13.6 eV, which have a large cross section. Thus, during the era 
of (ordinary) recombination, these photons have a short mean free path against absorption 
by neighboring ground-state hydrogen atoms. Thus, at this epoch the universe is optically 
thick to Lyman limit photons. Consequently, the overwhelming majority of recombinations 
to the ground state in one atom lead to a reionization of a neighboring atom, and there is 
no net change in the number of atoms. 

Ordinary recombination therefore proceeds via transitions initially to excited states, par- 
ticularly the n = 2 first excited state, and then to the ground state. However, the cosmic 
plasma is also optically thick to 2P — > IS Lyman-a photons, and so the transition to the 
ground state is dominated by the much slower two-photon 2S — > IS transition. Because of 
these effects, ordinary hydrogen recombination is not instantaneous, but delayed due to the 
'bottleneck' of the large optical depth for Lyman series photons. 

Our NLSP case of pX~ — » (pX) and 4 HeX _ — > ( 4 HeX) recombination is controlled by 
the same underlying atomic physics, and under conditions of a similar matter-to-photon 
ratio. It thus is worthwhile to check whether we expect similar effects. 

For the hydrogen and NLSP recombination, we are interested in the absorption of Lya 
and 'XLya:' photons, respectively, by ground-state atoms: 

71^2 + (Ap)i ->• (Ap) 2 , (19) 

where the 'anion' A G (e, X) corresponds to ordinary and NLSP recombination, respectively. 
Note, however, that in the NLSP case the proton, being lighter, plays the role of the electron 
in setting the relevant reduced mass p, so p = m e for hydrogen and p = m p for NLSP. 
On the other hand, the atomic mass m p + m A — B A is well approximated by m = m p for 
hydrogen and m = mx for NLSP. 

The XLya photon optical depth against absorption by (pX) atoms is 

T a (Ap) = cr a (Ap) n A d hor w a a (Ap) n A t , (20) 

so the ratio of hydrogen and NLSP optical depths is 

T a (Xp) <T a (Xp) Tlx t hhn 



T a (ep) o a (ep) n e t cmh 
The XLya resonance cross section is, in the notation of ref |123j . 



(21) 



where A a oc 1/E a ~ (a 2 /i) _1 , and r2 P ^i s oc p is the decay rate. Thermal broadening 
dominates over this width, with 5uj/lo ~ vt/c ~ yjT/m where m is the atomic mass. Thus 
we have an effective mean cross section 

r 2 r 2 
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where z is the redshift. The appropriate number densities are the physical, not comoving, 
values, and are set by Ti e ~ Tip-, and by n x = Y x n p , with n p ~ n B oc (1 + z) 3 . Finally, in the 
ordinary recombination case we have z cia b ~ 1000 and t cm b ~ 400, 000 yr ~ 10 13 sec, whereas 
in the NLSP case we have Zbbn ~ 4 x 10 8 and tbb 

n ~ 100 sec. 

Putting the above information together, we have 

r a (Xp) _ (m e \ 2 (m x \ y ( 1 + Zbbn \ 2 tbbn ^ 



TaifiP) J \ ?Tlp / \ ^ ^cmb / tcmb 

- 5 x 10- 7 ( *L) ( m * ) . (25) 
VO.Oiy VlOOGeV/ v ; 

Thus we see that the XLya optical depth at NLSP recombination is much smaller than 
that of ordinary recombination. However, the optical depth for ordinary recombination is 
enormous, r(ep) ~ 10 9 , so 

We find an optical depth against (pX) absorption that is much larger than unity. Hence it 
would seem that this effect could also be important for the NLSP recombination case. 

However, in the case of NLSP recombination there is an additional process to be con- 
sidered that has no analogue in ordinary recombination, namely the Compton scattering of 
XLya photons on free electrons and positrons. In ordinary recombination, electrons act as 
both the dominant photon scattering agent when they are unbound, and as the negatively- 
charged partners in the bound states. However, in our (AX~) case, these roles are now 
separated to and X~ respectively. The optical depth against electron scattering can be 
estimated using the ordinary Thomson cross section ct t : 

Te, = < b „>Ttbbn > <V T t bbn (27) 

^ 5 x io 7 I Tb ^ ) . (28) 

Viookevy V ; 

This is a lower bound, because we have used the net electron number n e)Iiet = n e - — n e + ~ n B , 
whereas pairs dominate the total e ± budget: 

n e - + n e + /m e \ 3 / 2 



1 ^ e ~m e/ T ^ 1 (2Q) 



n B 

down to T ~ m e /ln77 _1 ~ m e /25 ~ 20 keV. Note that during ordinary recombination, 
the optical depth against Thompson scattering drops below r e7 ~ 1, and cannot compete 
successfully with resonant Lya scattering. 

We see that an XLja photon will typically suffer at least ~ T ei /r a (Xp) ~ 500 Compton 
scatterings before encountering a bound state that it could reionize. Each of these scatter- 
ings degrades the photon energy, pulling it out of resonance. If the scattering were off a 
nonrelativistic electron, the photon would lose energy according to the Compton formula 

E - = TT^TT TV (30) 

1 H — 1 1 — cos 6* 
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and we would expect an approximate mean energy loss per scattering of 



^l^El. (31) 

Ery Tile 



The Lyman photons of interest have E(XLya) = 3/4 B (AX), and we have B( 4 B.eX) = 348 
keV and B(pX) = 25 keV; each species recombines at roughly T ~ B / lnr^ -1 ~ B/25. 

Thus 4 He recombination occurs when pairs are abundant, whereas protons recombine 
when pairs have completely annihilated. In either case, the Compton opacity dominates the 
resonance opacity, so that Lyman photons scatter many times before encountering a bound 
state. Moreover, in the first scattering the Lyman photons suffer energy losses AE/E ~ 0(1) 
for ( 4 HeX) and AE/E ~ (9(10 _1 ) for (pX). The Lyman photons are thus thermalized 
rapidly, long before they interact with any ground state atoms. We conclude that NLSP 
recombination to the ground state can occur unimpeded, unlike the case of ordinary recom- 
bination. 



B Specification of the Supersymmetric Model Frame- 
work 

In this paper we analyze the possible implications of bound states of massive metastable par- 
ticles on Big-Bang Nucleosynthesis in the context of the minimal supersymmetric extension 
of the Standard Model (MSSM). In this model, there is a supersymmetric partner for each 
Standard Model particle, and there are two Higgs doublet supermultiplets linked via a mix- 
ing parameter \x. The interactions are restricted to the same gauge and Yukawa interactions 
as in the Standard Model, so the quantity R = (— \) 3B + L + 2S [ s conserved multiplicatively, 
where B and L are the baryon and lepton numbers, respectively, and S is the spin. As a 
consequence of R conservation, the lightest supersymmetric particle (LSP) is stable, and a 
candidate for cosmological dark matter. 

We further assume the presence of soft supersymmetry-breaking fermion masses mu2, 
scalar masses mo and trilinear parameters A which are each universal at the grand unifi- 
cation scale, a framework known as the constrained MSSM (CMSSM) [95llll7j . In addition 
to m^jffio and A , we treat the ratio of Higgs vacuum expectation values, tan/3, as a 
free parameter. Motivated by the apparent discrepancy between the experimental mea- 
surement [118] of the anomalous magnetic moment of the muon, — 2, and theoretical 
calculations within the Standard Model [119] we assume that the MSSM Higgs mixing pa- 
rameter fj, is positive. As specified, the CMSSM includes no prediction for the mass of the 
gravitino, m 3 / 2 , which we treat as a free and independent parameter. 

We consider here the case in which the gravitino is the LSP, and the next-to-lightest 
supersymmetric particle (NLSP) is the spartner of one of the Standard Model particles. 
Possible candidates include the lightest neutralino \i t ne lighter stau slepton fj., the spartner 
of the right-handed electron or muon, e# or /2r, or the lighter stop squark t\. In the CMSSM 
as described above with a gravitino LSP, the most generic of these candidates for the NLSP 
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are the lightest neutralino \ and the lighter stau slepton fx, and the latter is the candidate 
we consider as an example of a charged metastable NLSP. 

C Three-Body Stau Decays 

We give here the matrix elements for f~ — > Gt~Z and f~ — > Gv T W~, which are the three- 
body stau decay processes that are most relevant for our study. In the Feynman diagrams 
for these two processes, the vertices involving the outgoing gravitino are given in Appendix 
A of |120j . and the vertices involving only the MSSM fields are given in |121] and |122j . In 
the context of the CMSSM, left-right mixing needs to be taken into account only for the 
third generation of sfermion fields. Neutrinos are treated the same way as in the Standard 
Model, i.e., as massless, purely left-handed neutrinos (and right-handed anti-neutrinos). 
Following Appendix B of [ID], we write the relation between the mass eigenstates, f 12 , and 
the interaction eigenstates, tl r, as 



r L \ _ / UfXL Uf2L \ ( f x 
TR J I U flR U f2 R ) 1 T 2 



(32) 



In order to be quite general, we give expressions for fj [j = 1,2), where the stau NLSP is 
the lighter of the two mass eigenstates. 

At tree-level, the Feynman diagrams contributing to fj~(px) —> G(j9 2 ) T~(p%) Z{p±) are the 
contact diagram and the t~ , f k (k = 1,2) and neutralino x° = 1 — 4) exchange diagrams. 
The partial matrix elements for each of these diagrams (suppressing spin and polarization 
indices) are 



^contact = lpn(p2) ( ~~iiT~ V~ ) [{ T f ~ Qf &w) U fjL P R + Q f SIB. 2 9 W Uf jR P L ] f] 



Mp cos 9w 

v(p 3 )e*M, (33) 



iM--^ T -*^GZT- 



i\/2~\ M-^+l + mr 



M P J (pi — P2 Y — m 2 - 

19 ^ l u [{Tf - Qf sin 2 Bw) Pr — Qf sin 2 9 W P L ] v(p 3 ) e* u (p 4 ), 

(34) 



cos 9 W 



iM 



1^2 



(u fkR p L -u fkL PR) {px-p±y 



(pi - Pi) - m 



cos 9 



w 



[{Tf - Q f sin 2 9 W ) U f3L U* fkL - Qf sin 2 9 w U fjR U* fkR ] 



(2px -Pi) u v(p 3 )el(p 4 ), 



(35) 
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+ (H R rr-\G R [f 4 ,'f}'f)F R 



% (ti-h + m xi) 



■ % 



(Pi -Ps) -m\ 
[ (C L U fjL + D L U fjR ) P L + (C R U fjR + D R U fjL ) P R ] v(p 3 ) e*M 



(36) 



where ip^{p2) represents the outgoing gravitino with momentum p 2 , v{ps) represents the 
outgoing chiral fermion with momentum p 3 , e* v {p^) is the polarization four- vector for the 
outgoing gauge boson with momentum p±, rf 11 is the flat-space Lorentz metric tensor, 
H L = m z (cosPN* k3 -sm(3N* k4 ), G L = cos6 w N* k2 - sm9 w N^ H R = H£, G R = G* L , 
C L = -(gm T N* k3 )/(V2m w cos/3), D L = V2gN* kl t&n9 w Q f , C R = C£, D R = ~V2g[N k2 Tf + 
tanOwNki [Qf — Ty) ] , T| = — 1/2, Qf — — 1 and N is the unitary matrix used to diagonal- 
ize the neutralino mass matrix (details can be found in |122j ). 

For fj{p\) — > G(p 2 ) v T (p?>) W~(p4), the Feynman diagrams contributing to this process, 
corresponding to eq. f )33|) . (IM]) . f )35|) and f )36|) . are the contact diagram, and the t~, v t and 
chargino \ k (k = 1, 2) exchange diagrams, respectively. The partial matrix elements are 
obtained by making the substitutions g — > gV2 cos 9w, Tf — > 1/2 and Qf — > in eq. ( 1331) - 
( )35|) . rrif k — > mc T , UfkR — > 0, Uf ul — > 1, U~ kR — > and U~ kL — > 1 in eq. (135]) since there is no 
right-handed sneutrino in the MSSM. Also, in eq. (1361) m^o — >■ m-±, and now the coefficients 

* * *k X-k 

are F L = V2m w cos /3C/£ 2 , G L = £7^, = \/2m w sin f3V k2 , = V*i, C L = D L = 0, 
Cr = (gm T Uk2) I (V2rn w cos /?) and _D R = —gllki, where U and V are the unitary matrices 
used to diagonalize the chargino mass matrix. 



References 

[1] D. Lindley, Astrophys. J. 294 (1985) 1. 

[2] J. R. Ellis, D. V. Nanopoulos and S. Sarkar, Nucl. Phys. B 259 (1985) 175. 

[3] D. Lindley, Phys. Lett. B 171 (1986) 235. 

[4] R. J. Scherrer and M. S. Turner, Astrophys. J. 331 (1988) 19. 

[5] M. H. Reno and D. Seckel, Phys. Rev. D 37 (1988) 3441. 

[6] S. Dimopoulos, R. Esmailzadeh, L. J. Hall and G. D. Starkman, Astrophys. J. 330, 545 
(1988); S. Dimopoulos, R. Esmailzadeh, L. J. Hall and G. D. Starkman, Nucl. Phys. B 
311 (1989) 699. 

[7] J. Ellis et al, Nucl. Phys. B 337, 399 (1992). 



39 



[9 
[10 

[11 
[12 

[13 
[14 

[15 

[16 

[it; 



[19 

[20 

[21 
[22 
[23 
[24; 
[25 



M. Kawasaki and T. Moroi, Prog. Theor. Phys. 93 (1995) 879 |arXiv:hep-ph/9403364] . 
M. Kawasaki and T. Moroi, Astrophys. J. 452, 506 (1995). 

E. Holtmann, M. Kawa saki, K. Kohri and T. Moroi, Phys. Rev. D 60, 023506 (1999) 
|arXiv:hep-ph/9805405] . 



K. Jedamzik, Phys. Rev. Lett. 84, 3248 (2000) |arXiv:astro-p h/9909445| . 

M. Kawasaki, K. K ohri and T. Moroi, Phys. Rev. D 63 (2001) 103502 
|arXiv:hep-ph/0012279] . 

K. Kohri, Phys. Rev. D 64 (2001) 043515 |arXiv:astro-ph/0103411] . 

R. H. Cyburt, J. R. El lis, B. D. Fields and K. A. Olive, Phys. Rev. D 67 (2003) 103521 
jarXiv: astro-ph/02 11258] . 

K. Jedamzik, Phys. Rev. D 70 (2004) 063524 [arXiv:astro-ph/04023~4"4 ]; K. Jedamzik, 
Phys. Rev. D 70 (2004) 083510 |arXiv:astro-ph/0405583l . ~ 

M. Kawasaki, K. K ohri and T. Moroi, Phys. Lett. B 625 (2005 ) 7 
|arXiv:astro-ph/0402490] ; Phys. Rev. D 71 (2005) 083502 |arXiv:astro-ph/0408426] . 

J. R. Ellis, K. A. O live and E. Vangioni, Phys. Lett. B 619, 30 (2005) 
|arXiv:astro-ph/0503023] . 

K. Kohri, T. Moroi and A. Yotsuyanagi, Phys. Rev. D 73, 123511 (2006) 
|arXiv:hep-ph/0507245l . 

D. G. Cerdeno, K. Y. Choi, K. Jedamzik, L. Roszkowski and R. Ruiz de Austri, JCAP 
0606, 005 (2006) larXiv:hep-ph/0509275| . 

K. Jedamzik, K. Y. Ch oi, L. Roszkowski and R. Ruiz de Austri, JCAP 0607, 007 (2006) 
|arXiv:hep-ph/0512044| . 

K. Jedamzik, Phys. Rev. D 74, 103509 (2006) [arXiv:hep-ph/0604251| . 



F. D. Steffen, JCAP 0609, 001 (2006) [arXiv:hep-ph/ 0605306|. 

M. Pospelov, Phys. Rev. Lett. 98 (2007) 231301 | arXiv:hep-ph/0605215] . 

K. Kohri and F. Takayama, Phys. Rev. D 76 (2007) 063507 | arXiv:hep-ph/0605243] . 

M. Kusakabe, T. Kajino and G. J. Mathews, Phys. Rev. D 74, 023526 (2006) 
|astro-ph/0605255]. 



[26] R. H. Cyb urt, J. R. Ellis, B. D. Fi elds, K. A. Olive and V. C. Spanos, JCAP 0611 
(2006) 014 |arXiv:astro-ph/0608562| . 



40 



[27] K. Hamaguchi, T. Hatsuda, M. Kamimura, Y. Kino and T. T. Yanagida, Phys. Lett. B 
650, 268 (2007) | hep-ph/0702274| [HEP-PH]]. 

[28] C. Bird, K. K oopmans and M. Pospelov, Phys. Rev. D 78, 083010 (2008) 
[hep-ph/0703096] . 

[29] M. Kawasaki, K. Kohri and T. Moroi, Phys. Lett. B 649, 436 (2007) |hep-ph/0703122) . 

[30] T. Jittoh, K. Kohri, M. Koike, J. Sato, T. Shimomura and M. Yamanaka, Phys. 
Rev. D 76, 125023 (2007) |arXiv:0704.29ll [hep-ph]]; Phys. Rev. D 78, 055007 (2008) 
[arXiv:0805.3389 [hep-ph]]; and T. Jittoh, K. Kohri, M. Koike, J. Sato, T. Shimomura 
and M. Yamanaka, Phys. Rev. D 82, 115030 (2010) |arXiv: 100 1.12171 [hep-ph]]. 



[31] K. Jedamzik, Phys. Rev. D 77, 063524 (2008) jarXiv:0707.2070l [astro-ph]]: and JCAP 
0803, 008 (2008) |axXiyi0im5153| [hep-ph]]. 

[32] D. Cumberbatch, K. Ichikawa, M. Kawasaki, K. Kohri, J. Silk and G. D. Starkman, 
Phys. Rev. D 76, 123005 (2007) |arXiv:0708.0095l [astro-ph]]. 

[33] J. Pradler and F. D. Steffen, Phys. Lett. B 666, 181 (2008) jarXiv:0710.22l3l [hep-ph]]. 

[34] M. Kusakabe, T. Kajino, R. N. Boyd, T. Yoshida and G. J. Mathews, Phys. Rev. 
D 76, 121302 (2007) jarXiv:0711.385l [astro-ph]]; and Astrophys. J. 680 (2008) 846 
jarXiv:0711.3858l [astro-ph]]. 

[35] M. Kawasaki, K. Kohri, T. Moroi and A. Yotsuyanagi, Phys. Rev. D 78, 065011 (2008) 
|arXiv:0804.3745l [hep-ph] ] . 

[36] M. Pospelov, J. Pradler and F. D. Steffen, JCAP 0811 (2008) 020 |arXiv:0807.4287l 
[hep-ph]]. 

[37] S. Bailly, K. Jedamzik and G. Moultaka, Phys. Rev. D 80 (2009) 063509 
|arXiv:0812.0788 l [hep-ph] ] . 

[38] M. Kusakabe, T. Kajino, T. Yoshida and G. J. Mathews, Phys. Rev. D 80, 103501 
(2009) |arXiv:0906.3516l [hep-ph]]. 

[39] K. Jedamzik and M. Pospelov, New J. Phys. 11, 105028 (2009) |arXiv:0906.2087l [hep- 
ph]]. 

[40] R. H. Cyburt, J. Ellis, B. D. Fields, F. Luo, K. A. Olive and V. C. Spanos, JCAP 0910 

(2009) 021 |arXiv:0907.5003l [astro-ph. CO]]. 

[41] M. Kusakabe, T. Kajino, T. Yoshida and G. J. Mathews, Phys. Rev. D 81, 083521 

(2010) jarXiv:1001.1410l [astro-ph. CO]]. 

[42] M. Pospelov and J. Pradler, Phys. Rev. D 82, 103514 (2010) jarXiv: 1006.41721 [hep-ph]] . 



41 



[43] R. H. Cyburt, J. Ellis, B. D. Fields, F. Luo, K. A. Olive and V. C. Spanos, JCAP 1010 
(2010) 032 [arXiv: 1007.41731 [astro-ph.CO]]. 

[44] M. Pospelov and J. Pradler, Phys. Rev. Lett. 106, 121305 (2011) [ arXiv:1010.4"079l 
[astro-ph.CO]]. 

[45] M. Pospelov and J. Pradler, Ann. Rev. Nucl. Part. Sci. 60, 539 (2010) jarXiv: 1011 .10541 
[hep-ph]]. 

[46] M. Kawasaki and M. Kusakabe, Phys. Rev. D 83, 055011 (2011) jarXiv: 1012.04351 [hep- 
ph]]. 

[47] D. A. Vasquez, A. Belikov, A. Coc, J. Silk and E. Vangioni, larXiv: 1208.04431 [astro- 
ph.CO]. 

[48] B. D. Fields, Ann. Rev. Nucl. Part. Sci. 61, 47 (2011) |arXiv:1203.3551 [astro-ph.CO]]. 

[49] G. Steigman, Ann. Rev. Nucl. Part. Sci. 57, 463 (2007) [ arXiv: 0712.11001 [astro-ph]]. 

[50] E. Komatsu et al. [WMAP Collaboration], Astrophys. J. Suppl. 192 (2011) 18 
[arXiv: 1001 .45381 [astro-ph.CO]]. 

[51] R. H. Cyburt, B . D. Fields and K. A. Olive, Astropart. Phys. 17, 87 (2002) 
[astro-ph/0105397] . 

[52] R. H. Cyburt, B . D. Fields and K. A. Olive, Phys. Lett. B 567, 227 (2003) 
[astro-ph/030243l] ; R. H. Cyburt, B. D. Fields and K. A. Olive, JCAP 0811, 012 
(2008) |arXiv:0808.28l8 [astro-ph]]. 

[53] A. Coc, E. Vangioni-Flam, P. Descouvemont, A. Adahchour and C. Angulo, Astrophys. 
J. 600, 544 (2004) [astro-ph/ 03 0"9480] . 

[54] R. H. Cyburt, B D. Fields and K. A. Olive, Phys. Rev. D 69, 123519 (2004) 
[astro-ph/031262 9|. 

[55] C. Angulo et al, Astrophys. J. 630, L105 (2005) |arXiv:astro-ph/0508454] . 

[56] R. N. Boyd, C. R. Brune, G. M. Fuller and C. J. Smith, Phys. Rev. D 82, 105005 (2010) 
jarXiv: 1008.08481 [astro-ph.CO]]. 

[57] R. H. Cyburt and M. Pospelov, Int. J. Mod. Phys. E 21 (2012) 1250004 |arXiv:0906.4373l 
[astro-ph.CO]]. 

[58] N. Chakraborty, B. D. Fields and K. A. Olive, Phys. Rev. D 83, 063006 (2011) 
|arXiv:1011.0722l [astro-ph.CO]]. 

[59] C. Broggini, L. Canton, G. Fiorentini and F. L. Villante, JCAP 1206, 030 (2012) 
[arXiv: 1202.52321 [astro-ph.CO]]. 



42 



[60] P. D. O'Malley, D. W. Bardayan, A. S. Adekola, S. Aim, K. Y. Chae, J. A. Cizewski, 
S. Graves and M. E. Howard et al, Phys. Rev. C 84, 042801 (2011). 



[61] O. S. Kirsebom and B. Davids, Phys. Rev. C 84, 058801 (2011) [ arXiv: 1109 .4690 [astro- 
ph.CO]]. 

[62] S. Vauclair,and C. Charbonnel, Ap. J. 502 (1998) 372 |arXiv:astro-ph/9802315l ; 
M. H. Pinsonneault, T. P. Walker, G. Steigman and V. K. Narayanan, Ap. J. 527 
(1998) 180 |[ arXiv:astro-ph/9803073] ; M. H. Pins onneault, G. Steigman , T. P. Walker, 
and V. K. Narayanan, Ap. J. 574 (2002) 398 | arXiv:astro-ph/0105439l; Q. Richard , 
G. Michaud and J. Richer, Astrophys. J . 619, 538 (2005) [arXiv:astro-ph/04 09672|; 
A. J. Korn et al, Nature 442, 657 (2006) |arXiv: astro-ph/ 060820 1| ; A. E. Garcia Perez, 
S., Inoue, W. Aoki, and S. G. Ryan, in Precision Spectroscopy in Astrophysics, Pro- 
ceedings of the ESO /Lisbon/ A veiro Conference, 9 (2008). 

[63] F. Spite, M. Spite, Astronomy & Astrophysics, 11 5 (1982) 357; S. G. Ryan , J. E. Norris 
and T. C. Beers, Astrophys. J. 523, 654 (1999) |arXiv:astro-ph/9903059] ; W. Aoki et 
al, Astrophys. J. 698, 1803 (2009) [arXiv: 0904. 14481 [astro-ph.SR]]: L. Sbordone et al, 
Astron. Astrophys. 522, 26 (2010) |arXiv:1003.45l0l [astro-ph. G A]]. 

[64] V.V. Smith, D.L. Lambert, and P.E. Nissen, Astrophys. J. 408, 262 (1993); Astrophys. 
J. 506, 405 (1998); L.M. Hobbs and J.A. Thorburn, Astrophys. J. Lett., 428, L25 
(1994); Astrophys. J. 491, 772 (1997); R. Cayrel, M. Spite, F. Spite, E. Vangioni- 
Flam, M. Casse, and J. Audouze, Astron. Astrophys. 343, 923 (1999); P. E. Nissen, 
M. Asplund, V. Hill and S. D'Odorico, Astron. Astrophys. 357, L49 (2000). 

[65] J. Melendez and I. Ramirez, Astrophys. J. 615, L33 (2004) [astro-ph/04093831 - 

[66] A. Hosford, A. E. G. Perez, R. Collet, S. G. Ryan, J. E. Norris and K. A. Olive, Astron. 
Astrophys. 493, 601 (2009) [ arXiv: 1004 .08631 [astro-ph.SR]]: A. Hosford, S. G. Ryan, 
A. E. G. Perez, J. E. Norris and K. A. Olive, Astron. Astrophys. 511, 47 (2010) 
|arXiv:0811.2"506 [astro-ph]]. 

[67] J. C. Howk, N. Lehner, B. D. Fields and G. J. Mathews, Nature 489, 121 (2012) 
|arXiv:1207.308H [astro-ph. CO]]. 

[68] C. Amsler et al [Particle Data Group], Phys. Lett. B667, 1 (2008). 

[69] I. Sick and D. Trautmann, Nucl. Phys. A 637 559 (1998). 

[70] D. H. Beck et al, Phys. Rev. C 30 1403 (1984); 

TUNL Nuclear Data Evaluation |http : //www . tunl . duke . edu/nu cldata/ 

[71] P. Mueller et al, Phys. Rev. Lett. 99 252501 (2007). 

[72] I. Sick, Phys. Rev. C 77 041302 (R) (2008). 

[73] R. Sanchez et al, Phys. Rev. Lett. 96 033002 (2006). 



43 



[74] W. Nortershauser et al, Phys. Rev. Lett. 102 062503 (2009). 

[75] M. Kamimura, Y. Kino and E. Hiyama, Int. J. Mod. Phys. A 24 (2009) 2076. 

[76] M. Pospelov. laFXiv:0712.0647 vl. (2007). 

[77] K. Arai, P. Descouvemont, D. Baye and W. N. Catford, Phys. Rev. C 68, 014310 (2003); 
K. Arai, Phys. Rev. C 69, 014309 (2004); also private communication cited in ref. [75]. 

[78] K. Varga and Y. Suzuki, Physical Review 52, 2882 (1995); K. Varga and Y. Suzuki, 
Comput. Phys. Commun. 106, 157 (1997). 

[79] S. Ali, and A. R. Bodmer, Nuclear Physics 80, 99 (1966). 

[80] S. Ali, and A. Afzal, II Nuovo Cimento 50, 355-358 (1967). 

11] D.R. Tilley, et al, Nucl. Phys. A636, 247 (1998). 

[82] G. R. Caughlan and W. A. Fowler, Atom. Data Nucl. Data Tabl. 40, 283 (1988). 

[83] C. Angulo, M. Arnould, M. Rayet, P. Descouvemont, D. Baye, C. Leclercq-Willain, 
A. Coc and S. Barhoumi et al, Nucl. Phys. A 656, 3 (1999). 

[84] R. A. Malaney and W. A. Fowler, Astrophys. J. 333, 14 (1988), Annals Phys. 192, 45 
(1989) 



[85] R. H. Cyburt, Phys. Rev. D 70, 023505 (2004) |astro-ph/04 01091|. 

[86] A. Coc, S. Goriely, Y. Xu, M. Saimpert and E. Vangioni, Astrophys. J. 744, 158 (2012) 
jarXiv: 1107.11171 [astro-ph.CO]]. 

]7] E. Matt, H. Pfander, H. Rieseberg, & V. Soergel, Physics Letters 9, 174 (1964). 

J. L. Feng, A. Rajar aman and F. Takayama, Phys. Rev. D 68 (2003) 063504 
|arXiv:hep-ph/0306024|; J. L. F eng, S. F. Su and F. Takayama, Phys. Rev. D 70 (2004) 
063514 [arXiv:hep-ph/0404198] . 

[89] J. L. Feng, S. Su and F. Takayama, Phys. Rev. D 70 (2004) 075019 
[arXiv:hep-ph/04042"31~l . 

[90] J. R. Ellis, K. A. Oli ve, Y. Santoso and V. C. Spanos, Phys. Lett. B 588 (2004) 7 
[arXiv:hep-ph/03122"62] . 

[91] J. R. Ellis, K. A. Olive, Y. Santoso and V. C. Spanos, Phys. Rev. D 70, 055005 (2004) 
[hep-ph/0405TT0| . 

[92] T. Sjostrand, P. Eden, C. Friberg, L. Lonnblad, G. Miu, S. Mrenna and E. Norrbin, 
Comput. Phys. Commun. 135 (2001) 238 |arXiv:hep-ph/0010017| . 



44 



[93] G. Aad et al. [ATLAS Collaboration], |M2£ivil207/[214| [hep-ex]; S. Chatrchyan et al. 
[CMS Collaboration]. [arXiv: 1207.72351 [hep-ex]]. 

[94] G. Degrassi, S. Heinemeyer, W. Hollik, P. Slavich and G. Weiglein, Eur. Phys. J. C 28 
(2003) 133 |arXiv:hep-ph/0212020]; S. Heinem eyer, W. Hollik and G. Weiglein, Eur. 
Phys. J. C 9 (1999) 343 larXiv:hep-ph/9812472]; S. Heinemeyer, W. H ollik and G. Wei- 
glein, Comput. Phys. Co mmun. 124 (2000) 76 [arXiv: hep-ph/ 9812320]; M. Fran k et al, 
JHEP 0702 (2007) 047 |arXiv:hep-ph/0611326] ; See |http : //www . f eynhiggs . de| 

[95] J. Ellis and K. A. Olive, Eur. Phys. J. C 72 (2012) 2005 |arXiv:1202.3"262l [hep-ph]]. 

[96] O. Buchmueller, R. Cavanaugh, A. De Roeck, M. J. Dolan, J. R. Ellis, H. Flacher, 
S. Heinemeyer and G. Isidori et al., arXiv: 11 12.3564 [hep-ph]; O. Buchmueller, R. Ca- 
vanaugh, M. Citron, A. De Roeck, M. J. Dolan, J. R. Ellis, H. Fl acher and S. Heinemeyer 
et al, larXiv: 1207.73151 [hep-ph]. 

[97] H. Baer, V. Barger and A. Mustafayev, Phys. Rev. D 85, 075010 (2012) [arXiv:1112.30T7l 
[hep-ph]]; A. Arbey, M. Battaglia, A. Djouadi, F. Mahmoudi and J. Quevillon, Phys. 
Lett. B 708, 162 (2012) [ arXiv:1112.3028l [hep-ph]]; L. Aparicio, D. G. Cerdeno and 
L. E. Ibanez, JHEP 1204, 126 (2012) |arXiv:1202.0"822l [hep-ph]]. 

[98] G. Aad et al. [ATLAS Collaboration], Phys. Lett. B 713, 387 (2012) [arXiv: 1204.07351 
[hep-ex]]; T. Aaltonen et al. [CDF Collaboration], Phys. Rev. Lett. 107, 239903 (2011) 
[Phys. Rev. Lett. 107, 191801 (2011)] [arXiv:1107.2"304l [hep-ex]]; updated results pre- 
sented at Aspen in Feb. 2012 by M. Rescigno, https://indico.cern.ch/getFile.py/ 
access?contribId=28&sessionId=7&resId =l&materialId=slides&conf Id=143360; 
S. Chatrchyan et al. [CMS Collaboration], Phys. Rev. Lett. 107, 191802 (2011) 
farXiv: 1107. 58341 [hep-ex] ] : R. Aaij et al. [LHCb Collaboration], Phys. Lett. B 699 (2011) 
330 |arXiv: 1103.21651 [hep-ex]]; Phys. Rev. Lett. 108, 231801 (2012) [arXiv:1203.4~493l 
[hep-ex]]; For an official combination of the ATLAS, CMS and LHCb results, see: 
ATLAS, CMS, and LHCb Collaborations, http : //cdsweb . cern. ch/record/1452186/| 
f iles/LHCb-CONF-2012-017 . pdf . 

[99] C. Bobeth, T. Ewerth , F. Kruger and J. Urban, Phys. Rev. D 64 (2001) 074014 
[arXiv:hep-ph/0104284|; A. Pedes, H . K. Dreiner and U. Nierste, Phys. Rev. Lett. 87 
(2001) 251804 larXiv:hep-ph/0108037|; C. Bobeth, A. J Buras, F. Kruger and J. Ur- 
ban, Nucl. Phys. B 630 (2002) 87 [arXiv:hep-ph/0 112305]; J. K. Mi zukoshi, X. Tata and 
Y. Wang, Phys. Rev. D 66 (2002) 115003 [arXiv:hep-ph/0208078l; C. Bobeth, T. Ew - 
erth, F. Kruger and J. Urban, Phys. Rev. D 66 (2002) 074021 [arXiv:hep-ph/0204225| ; 
A. J. Buras, P. H. Chankowski, J. Rosiek and L. Slawianowska, Nucl. Phys. B 659 
(2003) 3 [arXiv:hep-ph/0210145| ; A. J. Buras, M. V. Carlucci, S. Gori and G. Isidori, 
JHEP 1010 (2010) 009 [arXiv:1005.5'3T0l [hep-ph]]. 

[100] R. Arnowitt, B. Dutta, T. Kamon and M. Tanaka, Phys. Lett. B 538 (2002) 121 
larXiv:hep-ph/0203069l ; S. Baek, P. Ko and W. Y. Song, Phys. Rev. Lett. 89 (2002) 



45 



2718 01 |arXiv:hep-ph/ 0205259'|; C. S. Huang and X. H. Wu, Nucl. Phys. B 657 (2003) 
304 [arXiv:hep-ph/0212220[ ~S. Baek, Y. G. Kim and P. Ko, JHEP 0502 (2005) 067 
|arXiv:hep-p h/0406033 ] ; J. R. Ellis, K. A. Olive and V. C. Spanos, Phys. Lett. B 624, 
47 (2005) |hep-ph/0504196]; J. R. Ellis, K. A. Olive, Y. Santoso and V. C. Spanos, 
JHEP 0605 (2006) 063 [hep-ph/0603136] ; S. Heinemeyer, X. Miao, S. Su and G. Wei- 
glein, JHEP 0808, 087 (2008) [arXiv:0805.2"359l [hep-ph]]: B. Dutta, Y. Mimura and 
Y. Santoso, Phys. Lett. B 706, 188 (2011) [arXiv: 1107.3 020 [hep-ph]]; A. G. Akeroyd, 
F. Mahmoudi and D. M. Santos, JHEP 1112, 088 (2011) [arXiv: 1108.30181 [hep-ph]]: 
F. Mahmoudi, S. Neshatpour and J. Orloff. larXiv:1205.18 45 [hep-ph]. 

[101] K. Kohri and Y. Santoso, Phys. Rev. D 79, 043514 (2009) jarXiv:0811. 11191 [hep-ph]]. 

[102] K. A. Olive, P. Petitjean, E. Vangioni and J. Silk, |arXjyJ2m,5701] [astro-ph.CO]. 

[103] S. Buries and D. Tytler, Astrophys. J. 499, 699 (1998) |arXiv:astro-ph/9712108| ; 
S. Buries and D. Tytler, Astrophys. J. 507, 732 (1998) larXiv:astro-ph/9712109| ; 
J. M. O'Meara, D. Tytler, D. Kirkman, N. Suzuki, J. X. Prochaska, D. Lubin and 



A. M. Wolfe, Astrophys. J. 552, 718 (2001 ) |arXiv:astro-ph/0011179| ; M. Pettini and 
D. V. Bowen, Astrophys. J. 560, 41 (2001) |arXiv:astro-ph/0104474[ ; S. A. Levshakov, 
M. Dessauges-Zavadsky, S. D'Odorico and P. Molaro, Astrophys. J. 565, 696 (2002) 
[astro-ph/0105529|]; D. Kirkman, D. Tytler, N. Suzuki, J. M. O'Meara and D. Lubin, 



Astrophys. J. Suppl. 149, 1 (2003) larXiv:astro-ph70 302006| ; J. M. O'Meara, S. Buries, 
J. X. Prochaska, G. E. Prochter, R. A. Bernstein and K. M. Burgess, Astrophys. J. 649, 
L61 (2006) [arXiv:astro-ph/0608302l ; M. Pettini, B. J. Zych, M. T. Murphy, A. Lewis 
and C. C. Steidel, MNRAS 391, 1499 (2008) (arXiv: 0805 .05941 [astro-ph]] : R. Srianand, 
N. Gupta, P. Petitjean, P. Noterdaeme and C. Ledoux, MNRAS 405, 1888 (2010) 
[arXiv: 1002.46201 [astro-ph.CO]]; M. Fumagalli, J. M. O'Meara and J. X. Prochaska, 
Science 334 1245 (2011) jarXiv:1111.2'334l [astro-ph.CO]]. 

[104] G. S igl, K. Jedamzik, D. N. Schramm and V. S. Berezinsky, Phys. Rev. D 52 (1995) 
6682 [arXiv:astro-ph/9503094] . 

[105] J. Geiss, in Origin and Evolution of the Elements, eds. N. Prantzos, E. Vangioni-Flam, 
and M. Casse, (Cambridge University Press, Cambridge), 89 (1993). 

[106] K. A. Olive, and E. Skillman, New Astro nomy, 6, 119 (2001); K . A. Olive and 
E. D. Skillman, Astrophys. J. 617, 29 (2004) |arXiv:astro-"pTi704 05588 1 . 

[107] E. Aver, K. A. Olive and E. D. Skillman, JCAP 1103, 043 (2011) [ arXiv:1012.2"385l 
[astro-ph.CO]]. 

[108] Y. I. Izotov, T . X. Thuan and G. Stasinska, Astrophys. J. 662, 15 (2007) 
|astro-ph/0702072| [ASTRO-PH]]. 

[109] E. Aver, K. A. Olive and E. D. Skillman, JCAP 1204, 004 (2012) [ arXiv: 11 12.37131 
[astro-ph.CO]]. 



46 



[110] M. Asplund, D. L. Lambert, P. E. Nissen, F. Primas and V. V. Smith, Astrophys. J. 
644, 229 (2006) |arXiv:astro-ph/05 10636] . 

[Ill] R. Cayrel et al, Astron. Astrophys. 473, L37 (2007) jarXiv:0708.3"8l9 [astro-ph]]; 
M. Steffen, R. Cayrel, P. Bonifacio, H. G. Ludwig and E. Caffau, IAU Symposium, 265, 
23 (2010) |arXiv:0910.59T7 [astro-ph. SR]]. 

[112] G. Steigman, B. D. Fields, K. A. Olive, D. N. Schramm and T. P. Walker, Astrophys. 
J. 415, L35 (1993); B P. Fields and K. A. Olive, New Astronomy, 4, 255 (1999) 
|arXiv:astro-ph/9811183|; E. Vangioni-Flam, M. Casse, R. Cayrel, J. Audouze, M. Spite, 



and F. Spite, New Astronomy, 4, 245 (1999) |arXiv:astro-ph/9811327|. 

[113] S. G. Ryan, T. C. Beers, K. A. Olive, B. D. Fields, and J. E. Norris Astrophys. J. Lett. 
530 (2000) L57 |arXiv: astro-ph /99052 1 1] . 

[114] J. I. G. Hernandez et al, Astron. Astrophys. 505, L13 (2009) |arXiv: 0909 .09831 [astro- 
ph.GA]]; L. Pasquini and P. Molaro, Astron. Astrophys. 307, 761 (1996); F. Thevenin, 
C. Charbonnel, J. A. d. Pacheco, T. P. Idiart, G. Jasniewicz, P. de Laverny and B. Plez, 
Astron. Astrophys. 373, 905 (20 01) [arXiv:astro-ph/0105166| ; P. Bonifacio et al, As- 
tron. Astrophys., 390, 91 (2002) [arXiv:astro-ph/0204332] ;P Bonifacio, Astron. Astro- 
phys. 395, 515 (2002) |arXiv:astro-ph/0209434] ; K. Lind, F. Primas, C. Charbonnel, 
F. Grundahl and M. Asplund, Astron. Astrophys. 503, 545 (2009) |arXiv:0906.2876l 
[astro-ph. SR]]. 

[115] A. M. Boesgaard, J. A. Rich, E. M. Levesque and B. P. Bowler, Astrophys. J. 743, 
140 (2011) jarXiv:1110.2823l [astro-ph. SR]]. 

[116] H. Ito, W. Aoki, S. Honda and T. C. Beers, Astrophys. J. 698, L37 (2009) 
|arXiv:0905.0950l [astro-ph. SR]]. 

[117] M. Drees and M. M. Nojiri, Phys. Rev. D 47 (1993) 376 arXiv:hep-ph/9207234]; 



G. L. Kan e, C. F. Kolda, L Roszkowski and J. D. Wells, Phys. Rev. D 49 (1994) 
6173 larXiv:hep-ph/9312~272] ; H. Baer and M. Brhlik, Ph ys. Rev. D 53 (19 96) 
597 |arXiv:hep-ph/9508321] ; Phys. Rev. D 57 (1998) 567 |arXiv:hep-ph/9706509] ; 
J. R. Ellis, T. Falk, K . A. Olive and M. Schmitt, Phys. Lett. B 388 (1996) 97 
|arXiv:hep-ph/9607292 |; Phys. Lett. B 413 (1997) 355 larXiv:hep-ph/9705444] ; J. R. El- 
lis, T. Falk, G. Ganis, K. A . Olive and M. Schmitt, Phys. Rev. D 58 (1998) 
09500 2 |arXiv:hep-ph/9801445] ; V. D. Barger and C. Kao, Phys. Rev. D 57 (1998) 
3131 [arXiv:hep-ph/9704403]; J. R. Ellis, T. Fa lk, G. Ganis and K. A. Olive, Phys. 
Rev. D 62 (2000) 075010 |arXiv:hep-ph/00041~69] ; H. Baer, M. Brhlik, M. A. Diaz, 
J. Ferrandis, P. Mercadante, P. Quintana and X. Tata, Phys. Rev. D 63 (2001) 
015007 larXiv:hep-ph/000502 7| ; J. R. Ellis, T . Falk, G. Ganis, K A. Olive and 
M. Srednicki, Phys. Lett. B 510 (2001) 236 |arXiv:hep-ph/01020~9~8] ; V. D. Barger 
and C. Kao, Phys. Lett. B 518 (2001) 117 |arXiv:hep-ph/0106189]; L. Roszkowski , 
R. Ruiz de Austri and T. Nihei, JHEP 0108 (2001) 024 [arXiv:hep-ph/ 01 06334] ; 



47 



A. Djouadi, M. Drees and J. L. Kneur, JHEP 0108 (2001) 055 |arXiv:hep-ph/0107316] ; 
U. Chattopadhyay, A . Corsetti and P. Nath, Phys. Rev. D 66 (2002) 035003 
|arXiv:hep-p h/020 10"0~l~]; J. R. Ell is, K. A. Olive and Y. Santoso, New Jour. Phys. 4 
(2002) 32 larXiv:hep-ph/0202110' l; H. Baer, C. Bal azs, A. Belyaev, J. K. Mizukoshi, 
X. Tata and Y. Wang, JHEP 0207 (2002) 050 |arXiv:hep-ph70205325 1 ; R. Arnowitt and 



B. Dutta, arXiv:hep-ph/0211417 J. R. Ellis, K. A. Olive, Y. Santoso and V. C. Spanos, 



Phys. Lett. B 565 (2003) 176 |arXiv:hep-ph/0303043|; H. Baer and C. Balazs, JCAP 
0305, 006 (2003) |arXiv:hep-ph/0303114]; A. B. Lahanas and D. V. Nanopoulos, Phys. 
Lett. B 568, 55 (2003) [arXiv:hep-ph/0303130]; U. Chattopadhya y, A. Corsetti and 
P. Nath, Phys. Rev. D 68, 0350 05 (2003) [arXiv:hep-ph/030320T] ; C. Munoz, Int. J. 
Mod. P hys. A 19, 3093 (2004 ) [arXiv:hep-ph"70 309346 1 ; R. Arnowitt, B. Dutta and 
B. Hu, arXiv:hep-ph/0310103] J. Ellis and K. A. Olive, in Particle Dark Matter, ed. 
G. Bertone, ppl42-163 jarXiv:1001.365T1 [astro-ph.CO]]. 

[118] [The Muon g-2 Collaboration], Phys. Rev. Lett. 92 (2004) 161802, 
[arXiv:hep-ex/0401008 1 ; G. Be nnett et al. [T he Muon g-2 Collaboration], Phys. 
Rev. D 73 (2006) 072003 [arXiv:hep-ex/0602035] . 

[119] D. Stockinger, J. Phys. G 34 (2007) R45 larXiv:hep-ph/0609168]; J. Miller , 
E. de Rafael and B. Roberts, Rept. Prog. Phys. 70 (2007) 795 [arXiv:hep-ph/0703049] ; 
J. Prades, E. de Rafael and A. Vainshtein, larXiv:0 901.0306 [hep-ph]; F. Jegerlehner 
and A. Nyffeler, Phys. Rept. 477, 1 (2009) jarXiv:0902.3360l [hep-ph]]; M. Davier, 
A. Hoecker, B. Malaescu, C. Z. Yuan and Z. Zhang, Eur. Phys. J. C 66, 1 
(2010) |arXiv:0908.4300l [hep-ph]]; J. Prades, Acta Phys. Polon. Supp. 3, 75 (2010) 
Jar2Gvi0909.254l [hep-ph]] ; T. Teubner, K. Hagiwara, R. Liao, A. D. Martin and D. No- 
mura, larXiv:1001.5 401 [hep-ph]; M. Davier, A. Hoecker, B. Malaescu and Z. Zhang, 
larXiv: 1010.41801 [hep-ph]. 

[120] F. Luo, K. A. Olive and M. Peloso, JHEP 1010, 024 (2010) jarXiv:1006.5"570l [hep-ph]]. 

[121] H. E. Haber and G. L. Kane, Phys. Rept. 117, 75 (1985). 

[122] J. F. Gunion and H. E. Haber, Nucl. Phys. B 272, 1 (1986) [Erratum-ibid. B 402, 567 
(1993)]. 

[123] P. J. E. Peebles, Principles of Physical Cosmology (1993), see equations (6.100)-(6.104) 
and the surrounding discussion. 



48 



